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Preface 
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Administration  (NASA)  Atmospheric  Effects  of  Aviation  Project  (AEAP),  and  managed 
through  the  sub-element  Emission  Characterization  at  NASA  Lewis  Research  Center 
(LeRC)  by  Mr.  Richard  Niedzwiecki,  NASA,  and  Dr.  Chowen  Chou  Wey,  Army 
Research  Laboratory  (ARL).  The  engine  test  was  conducted  at  NASA  LeRC,  Cleveland, 
OH.  It  was  accomplished  by  the  team  efforts  of  many  engineers  and  scientists  from 
government  and  private  organizations  including  NASA  LeRC,  Dynacs  Engineering 
Company,  Inc.  (Dynacs),  Pratt  &  Whitney  (P&W),  Aerodyne  Research  Inc.  (ARI),  Air 
Force  Research  Laboratory  (AFRL),  and  University  of  Missouri-Rolla  (UMR). 


The  main  objective  of  this  report  is  to  document  the  results  of  a  set  of  parametric  studies  of 
the  effect  of  the  fuel  sulfur  content  on  particulate  emissions.  The  primary  authors  were: 


Dr.  John  O.  Ballenthin 
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Dr.  Thomas  M.  Miller 
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Dr.  Philip  D.  Whitefield 
Dr.  Joda  C.  Wormhoudt 


AFRL 
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UMR 

ARI 
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AFRL 

Dynacs/NASA  LeRC 
ARL/NASA  LeRC 
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ARI 


This  task  is  an  accumulation  of  contributions  from  individuals  acknowledged  below  and 
other  persons  too  numerous  to  list.  Although  participants  worked  as  a  team  with  contributions 
beyond  individual  responsibilities,  primary  functions  associated  with  each  were: 


Project  Management 
Principal  Investigator 


Mr.  Richard  W.  Niedzwiecki 
Dr.  Chowen  Chou  Wey 
Dr.  Chowen  Chou  Wey 
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ARL/NASA  LeRC 
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1.  INTRODUCTION 


The  overall  objective  of  the  NASA  Atmospheric  Effects  of  Aviation  Project  (AEAP)  is  to 
develop  scientific  bases  for  assessing  atmospheric  impacts  of  the  exhaust  emissions  by  both 
current  and  future  fleets  of  subsonic  and  supersonic  aircraft.  Among  the  six  primary  elements  of 
the  AEAP  is  Emissions  Characterization.  The  objective  of  the  Emission  Characterization  effort  is 
to  determine  the  exhaust  emission  constituents  and  concentrations  at  the  engine  exit  plane. 

In  light  of  the  recently  obtained  in-flight  measurements  of  the  composition  of  the  plume  of  a 
Concorde  aircraft,  as  described  in  a  paper  by  Fahey  et  al.  (1995)  and  in  NASA’s  AES  A  report 
(Stolarski  et  al.,  1995),  the  potentially  important  role  of  aerosol  and  aerosol  precursor  emissions 
is  recognized.  The  particulate  measurements  obtained  in  these  flight  measurements  indicate  a 
high  degree  of  sulfur  dioxide  (SO2)  oxidation  to  condensed  sulfate.  The  data  also  suggest  that 
SO2  oxidation  by  hydroxyl  radicals  (OH)  is  not  the  dominant  cause  of  the  conversion. 

Particulates  directly  emitted  by  jet  aircraft  are  mostly  soot  with  traces  of  heavy  unburned 
hydrocarbons.  Volatile  aerosols  are  thought  to  be  formed  in  the  exhaust  of  the  engines  as  a  result 
of  emissions  related  to  the  sulfur  impurity  in  aviation  fuel.  Sulfur-containing  aerosols  are  directly 
involved  in  the  partitioning  of  reactive  species  controlling  the  abundance  of  ozone  in  the 
stratosphere  where  supersonic  aircraft  principally  operate.  This  impact  is  maximized  if  the 
sulfate  is  formed  within  the  exhaust  plume.  In  the  troposphere,  subsonic  emissions  can  cause  the 
formation  of  contrails  that,  in  turn,  can  also  cause  changes  in  background  aerosol  properties 
which  may  affect  cloud  cover  and  cloud  optical  properties.  Changes  in  ozone  and  cloud 
properties  are  both  elements  of  climate  change  and  hence  can  influence  the  environment  on  long 
time  scales. 

Although  the  importance  of  aerosols  and  their  precursors  is  now  well  recognized,  the 
characterization  of  current  subsonic  engines  for  these  emissions  is  far  from  complete. 
Furthermore,  since  the  relationship  of  engine  operating  parameters  to  aerosol  emissions  is  not 
known,  extrapolation  to  untested  and  unbuilt  engines  necessarily  remains  highly  uncertain.  This 
engine  test,  as  well  as  a  parallel  flight  measurement,  attempts  to  address  both  issues  by 
expanding  measurements  of  aerosols  and  aerosol  precursors  with  fuels  containing  different  levels 
of  fuel  sulfur  content.  Aviation  fuels  with  a  range  of  ftiel  sulfur  were  procured  for  both  sampling 
venues. 

The  specific  objective  of  this  engine  test  is  to  obtain  a  database  of  sulfur  oxides  emissions  as 
well  as  the  non-volatile  particulate  emission  properties  as  a  function  of  fuel  sulfur  and  engine 
operating  conditions.  The  database  of  the  non-volatile  particulates  emission  properties  is  to  be 
used  as  a  comparative  baseline  with  subsequent  flight  measurement.  A  database  of  volatile  and 
non-volatile  particulates  and  sulfur  oxides  is  needed  to  validate  engine  and  exhaust  plume 
models. 

2.  TEST  FACILITY 

The  engine  used  in  this  test  was  a  Pratt  &  Whitney  F100-200E  turbofan  engine.  The  normally 
variable  convergent-divergent  type  exhaust  nozzle  was  mechanically  locked  into  position  so  that 
a  throat  area  of  0.2787  m2  (3.0  ft2)  would  be  maintained  throughout  the  entire  envelope  of  test 
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conditions.  The  fan  bypass  air  is  combined  with  the  engine  core  flow  upstream  of  the  exhaust 
nozzle,  with  no  intentional  mixing  of  the  two  streams  of  air  being  performed  before  the  gases 
exit  the  nozzle.  The  engine  is  representative  of  commercial  type  engines  in  the  89  -  133  KN 
(20,000  -  30,000  lbf)  thrust  class.  It  has  an  engine  pressure  ratio  of  25  with  a  bypass  ratio  of 
0.75.  Test  conditions  ranged  from  idle  to  military  power  settings  at  altitudes  from  sea  level  up  to 
16,764m  (55,000fit).  The  augmentor  was  locked  out  during  the  entire  test.  The  simulated  flight 
Mach  number  was  0.0  at  sea  level  and  0.8  for  all  other  altitude  conditions.  Aviation  fuels  that 
met  or  exceeded  standard  specification  were  used  throughout  the  test. 

As  shown  in  Figure  1,  the  engine  was  installed  into  the  NASA  LeRC  PSL-4  test  cell  on  a  multi¬ 
directional  thrust  measuring  stand.  Air  is  supplied  to  the  engine  through  0.8826  m  (34.75  in) 
diameter  inlet  ducts  which  are  directly  connected  to  the  front  of  the  engine.  Before  the  air  enters 
the  engine,  the  pressure  and  temperature  are  set  to  match  the  conditions  that  the  engine  would 
see  at  the  desired  altitude  and  flight  Mach  number.  Downstream  of  the  engine  exhaust  nozzle  is 
the  gas  sampling  rake  followed  by  the  1.83m  (72  in)  diameter,  water  cooled  exhaust  collector.  A 
circular  cooling  air  torus  is  located  in  the  forward  portion  of  the  test  cell  and  provided  9.09  to 
13.64  kg/sec  (20  to  30  lb/sec)  of  cooling  air  maintaining  a  temperature  below  37.8  °C  ( 100  °F). 

In  order  to  protect  the  engine  and  test  cell,  the  sea  level  static  test  points  were  run  at  total 
temperatures  and  total  pressures  equal  to  a  normal  day  at  sea  level,  but  the  altitude  tank  was  at  a 
pressure  equivalent  to  an  altitude  of  3,048  m  (10,000  ft).  These  types  of  conditions  will  yield 
engine  performance  very  similar  to  those  at  a  true  static  sea  level  condition. 

The  system  used  on  this  program  to  provide  a  hot,  pressurized  exhaust  gas  sample  to  the  gas  and 
particle  analyzers  consisted  of  a  cruciform  type  rake  located  behind  the  engine  exhaust  nozzle, 
which  held  the  sampling  probes,  a  distribution  box  and  a  pump  box.  The  vertical  arms  of  the 
cruciform  rake  contain  only  gas  sampling  probes.  The  horizontal  arms  contain  alternating  Mach 
Number/Flow  Angularity  (MFA),  stagnation  temperature,  and  gas  sampling  probes.  The  gas 
sampling  probes  were  designed  to  meet  emissions  certification  test  requirements.  They  were 
spaced  equally  across  the  arm,  and  could  be  ganged  (valved  to  a  common  line)  or  sampled 
individually.  Individual  probe  sampling  allows  spatial  mapping  of  engine  exit  plane  core,  bypass, 
and  intermediate  mixed  flow  regimes. 

The  rake  system  was  mounted  on  a  support  structure  positioned  between  the  engine  exhaust 
nozzle  and  the  inlet  to  the  exhaust  collector.  The  tips  of  the  sampling  probes  were  positioned 
0.2144  m  downstream  of  the  exhaust  nozzle  exit  plane.  Details  of  the  sampling  rake  along  with 
the  sampling  probe  design  is  described  in  Appendix  A. 

3.  TEST  PARAMETERS 

Aviation  fuel  (Jet  A)  with  a  range  of  fuel  sulfur  was  procured.  Low  and  high  sulfur  values  are 
limited  by  commercially  available  fuels  and  by  fuel  specification  limits  of  0.3%  by  weight.  Low- 
sulfur  fuel  used  had  a  sulfur  content  of  approximately  18  parts  per  million  (ppm),  and  the  high- 
sulfur  fuel  of  1 1 13  ppm.  An  intermediate  level,  152  ppm,  was  set  by  mixing  the  low  and  the  high 
sulfur  fuels.  A  military  fuel  JP-8  +  100  with  336  ppm  sulfur  was  also  tested.  Fuel  samples  were 
sent  to  both  Wright  Patterson  AFB  and  a  private  laboratory  to  be  comprehensively  analyzed  for 
each  batch.  Details  of  these  analyses  are  listed  in  Tables  1  and  2.  Sample  #1  is  the  low  sulfur  fuel 
with  corrosion  inhibitor.  Sample  #2  is  the  low  sulfur  fuel  with  corrosion  inhibitor,  conductivity 
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additive,  and  icing  inhibitor.  Sample  #3  is  the  high  sulfur  fuel  with  corrosion  inhibitor.  Sample 
#4  is  the  mixed  medium  sulfur  fuel  with  corrosion  inhibitor,  conductivity  additive,  and  icing 
inhibitor.  Sample  #5  is  the  JP8+100.  Fuel  additives  used  in  this  test  series  are  listed  in  Table  3. 

For  the  first  2  nights  of  the  test  program,  the  gas  samples  were  drawn  from  individual  sample 
probes  in  order  to  determine  the  distribution  of  the  exhaust  emissions  across  the  engine  exit 
plane.  This  also  provided  an  opportunity  to  calibrate  all  the  emission  measurement  systems. 
Individual  samples  were  drawn  from  the  vertical  rake  the  first  night,  then  the  horizontal  rake  the 
second  night.  Figures  2  and  3  show  the  C02  measured  from  individual  probes  as  a  typical 
emission  distribution. 

Based  on  these  distributions,  the  3  inner  most  sampling  probes  from  the  engine  centerline  of  each 
arm  were  ganged  together  for  each  measurement  system  for  the  third  and  subsequent  test  nights. 
Both  of  the  horizontal  arms  are  ganged  together  for  the  LeRC  system.  The  bottom  vertical  rake 
was  dedicated  to  the  UMR/MASS  system.  The  top  vertical  rake  was  dedicated  to  the 
AFRL/CIMS  system. 

Table  4  listed  the  test  matrix  which  was  set  by  parametrically  varying  the  combustor  inlet 
temperature  (T3)  between  idle  and  maximum  power  setting  at  simulated  SLS  and  up  to  five  other 
altitudes  for  each  fuel.  Figure  4  graphically  illustrates  the  test  conditions  as  a  function  of 
combustor  pressure  and  temperature. 

4.  EMISSION  MEASUREMENT  SYSTEMS 

Four  diagnostic  systems,  extractive  and  non-intrusive  (optical),  were  assembled  for  the  gaseous 
and  particulate  emissions  characterization  measurements  study.  These  systems  are  briefly 
described  below.  Details  of  the  respective  techniques  and  database  established  from  this  test  are 
described  in  the  appendices. 

4.1  NASA  Extractive  Gaseous  Emission  Measurement  System  and  Smoke  Meter 


The  extractive  gaseous  emissions  measurement  system  contains  an  array  of  analyzers 
dedicated  to  examining  the  concentrations  of  specific  gases  and  a  standard  SAE  smoke 
meter.  This  system  is  essentially  identical  to  systems  used  for  commercial  aircraft  engine 
emissions  certification  testing,  except  that  probes  are  equally  spaced  in  the  present  study, 
affecting  the  spatial  averages  of  the  exhaust.  The  gaseous  emissions  system  includes 
analyzers  for  measurement  of  CO,  C02,  NO,  NOx,  02,  total  unburnt  hydrocarbons  (THC), 
and  S02.  Details  of  each  analyzer  are  described  in  Appendix  B. 

Spatial  profiles  were  obtained  by  sampling  individual  probes  of  each  arm  of  the  rake 
system.  Comparisons  of  vertical  and  horizontal  species  concentration  and  temperature 
profiles  from  the  first  two  days  of  test  series  confirmed  symmetry  and  gave  confidence  in 
representative  sampling  from  either  arm.  Three  inner-most  probes  from  both  horizontal 
arms  were  ganged  for  this  system. 
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4.2  UMR  Extractive  Particulate  Measurement  System 


Appendix  C  contains  details  on  the  particulates  and  aerosols  characterization  system 
employed  in  this  study  -  the  University  of  Missouri-Rolla  Mobile  Aerosol  Sampling 
System  (MASS).  A  brief  overview  of  the  system  is  given  here.  For  this  test  the  UMR 
trailer-based  MASS  facility  was  employed.  The  trailer  was  located  immediately  adjacent 
to  the  NASA  LERC  PSL  building  and  was  connected  to  the  altitude  chamber  through 
30m  long  heated  sample  and  diluent  lines.  Particulate  concentrations  from  single  and 
ganged  (multiple)  probe  sources  and  were  monitored  in  real  time,  with  typically  a  1  Hz 
sampling  frequency,  using  commercially  available  condensation  nucleus  counters 
(CNC's)  (TSI  and  Met-One  models:  3025  and  1 105-5,  respectively).  Size  distributions 
were  determined  both  in  real  time  and  from  sample  tanks  for  particulate  diameters 
ranging  from  7  -  250  nm  using  differential  mobility  analysis.  Hydration  properties  (e.g. 
soluble  mass  fraction)  were  derived  from  tank  samples  using  a  tandem  differential 
mobility  analysis  system.  A  needle  to  grid  electrostatic  precipitator  was  used  to  collect 
particulate  samples  on  electron  microscope  grids.  Comprehensive  data  presented  in 
Appendix  D  are  summarized  in  section  5.3. 

4.3  AFRL  Extractive  Sulfur  Oxides  Measurement  System 

A  chemical  ionization  mass  spectrometer  (CIMS)  from  the  Air  Force  Research 
Laboratory  at  Hanscom  AFB  was  used  to  measure  levels  of  SCL  and  HNO3  in  the  F-100 
jet  engine  exhaust.  Exhaust  gases  were  sampled  through  three  extractor  probes  on  the 
upper  part  of  the  vertical  rake,  26  cm  beyond  the  exit  plane  of  the  engine.  The  inlet  lines 
were  joined  at  a  point  1  m  from  the  probes,  where  the  sampled  gases  entered  an  0.64-cm 
OD,  7.3-m  long  stainless  steel  line  which  transported  the  exhaust  effluent  to  the  flow  tube 
of  the  CIMS.  The  sampling  line  was  heated  to  150  °C  to  prevent  condensables  from 
adhering  to  the  walls  of  the  sampling  line.  The  stainless  steel  lines  (including  the  three 
probe  lines)  were  coated  with  siloxane-covered  glass  to  minimize  sticking  of  HNO3  to  the 
walls.  The  exhaust  effluent  entered  the  CIMS  flow  tube  through  a  0.32-cm  OD,  7.5-cm 
long  capillary,  which  restricted  the  effluent  flow  to  4-15  slm  (standard  liters  per  minute). 
The  flow  tube  utilized  a  fast  flow  (120  slm)  of  N2  gas  at  37  Torr  pressure,  diluting  the 
engine  effluent  by  typically  a  factor  of  6. 

Detection  of  the  trace  gases  in  the  engine  effluent  was  made  by  reacting  the  diluted 
engine  effluent  with  CO3'  ions  in  the  flow  tube  over  a  path  length  of  2.5  cm.  The 
chemistry  involved  in  detecting  SO2  and  HNO3  is: 

CO3  +  SO2  — >  SO3  +  CO2, 
followed  in  less  than  a  microsecond  by: 

SO3  +  O2  +  M  — >  SO5  +  M, 

where  M  is  a  third  body  (Ni  or  O2).  Thus,  the  appearance  of  ion  signal  at  1 12  amu  (SO.O 
signifies  the  presence  of  SO2  in  the  engine  exhaust.  Similarly, 

CO3  +  HNO3  +  M  — >  CO3XHNO3)  +  M, 
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and  thus,  appearance  of  ion  signal  at  123  amu  [CO3XHNO3)]  signifies  the  presence  of 
HNO3  in  the  engine  exhaust.  Known  flow  rates  of  calibration  gases  (SO2  and  HNO3) 
were  periodically  added  during  each  engine  test  condition  (altitude,  power  level,  and  fuel 
type).  Tests  were  carried  out  to  determine  loss  rates  in  the  sampling  line,  by  injecting 
calibrant  gases  consecutively  at  four  points  along  the  line  (see  Appendix  E).  The 
combination  of  calibrant  ion  signal,  calibrant  flow  rate,  engine  exhaust  ion  signal,  and 
total  engine  effluent  flow  rate,  give  the  mixing  ratio  (fractional  concentration  of  SO2  or 
HNO3)  in  the  engine  exhaust.  Emission  indices  (grams  of  SO2  or  NO2  equivalent  per 
kilogram  of  fuel,  respectively)  were  evaluated  from  the  respective  mixing  ratios  using 
NASA’s  data  for  the  percentage  of  CO2  combustion  product  in  the  exhaust.  Details  of  the 
method,  instrument,  and  data  analysis  are  given  in  Appendix  E. 

4.4  ARI  Non-Intrusive  Gaseous  Emissions  Measurement  System 

Aerodyne  Research,  Inc.  (ARI)  used  infrared  tunable  diode  laser  (IR-TDL)  absorption  to 
measure  S02,  SO3,  NO,  H2O  and  C02  along  an  optical  path  through  the  center  of  the 
plume.  Two  lasers  were  operated  simultaneously,  with  the  two  spatially  coincident 
beams  focused  on  the  same  detector.  Multiplex  detection  was  achieved  by  turning  off 
each  laser  while  the  other  laser  was  being  spectrally  scanned.  The  data  acquisition  rate 
resulted  in  accumulation  of  spectra  from  both  lasers  at  a  600  Hz  rate.  The  laser  beams 
entered  and  exited  the  test  cell  through  a  calcium  fluoride  window  with  a  1-deg  wedge  to 
suppress  interference  fringes.  The  optical  system  inside  the  test  cell  allowed  for  either  14 
passes  through  the  exhaust,  using  multi-pass  mirrors  inside  purged  protective  boxes,  or 
two  passes,  using  a  retro-reflector.  Some  readjustment  of  the  optical  system  during  the 
test  was  required,  due  to  motion  of  the  test  cell  with  changes  in  simulated  altitude. 
During  early  test  days,  spectra  were  taken  of  several  NO  and  CO2  lines,  in  both  two-pass 
and  multi-pass  configurations.  For  some  low  sulfur  cases,  and  for  all  observations  of 
fuels  with  higher  sulfur  loadings,  the  two  spectral  regions  observed  contained  lines  of 
S02  and  SO3,  respectively.  Both  regions  also  contained  H2O  lines,  allowing  direct 
calculation  of  sulfur  oxide  emission  indices  using  the  known  water  emission  index.  The 
spectral  data,  analysis  procedures,  and  results  are  presented  in  Appendix  F. 

5.  SUMMARY  OF  RESULTS 

5.1.  Gaseous  Emission  Data 


Multiple  data  points  acquired  during  steady- state  conditions  were  relatively  constant  over 
time  and  were  therefore  averaged  and  reported  as  a  single  value.  Dates  and  clock  times 
are  included  in  the  table  to  document  approximate  duration  per  each  steady-state  test 
condition,  and  to  correlate  with  data  from  other  measurement  systems.  The  full  set  of 
data,  including  engine  parameters,  smoke  number  and  gaseous  emissions,  is  listed  and 
discussed  in  Appendix  C. 
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5.2.  Comparison  of  Extractive  and  Non-Intrusive  Sulfur  Oxides  Data 

SO2  emissions  were  measured  by  three  different  techniques:  extractive  sampling  using 
UV  absorption  (LeRC)  and  CIMS  (AFRL)  as  well  as  non-intrusive  in  situ  infrared  TDL 
absorption  (ARI).  The  SO2  emissions  indices  for  all  three  techniques  agree  within  their 
respective  uncertainty  estimates.  Emission  indices  averaged  over  all  T3  values  are  plotted 
in  Figure  5  as  a  function  of  altitude  for  the  three  highest  fuel  sulfur  loadings  and  all  three 
measurement  techniques.  Table  5  gives  emission  indices  averaged  over  all  altitudes  and 
power  settings.  Figure  5  shows  that  within  the  error  limits,  there  is  no  trend  with  altitude. 
A  similar  lack  of  any  discernible  trend  applies  to  plots  as  a  function  of  T3.  Consideration 
of  any  column  in  Table  5  might  suggest  a  trend  with  sulfur  loading  for  a  particular 
measurement  technique,  but  a  comparison  among  all  three  techniques  indicates  that  no 
trend  exists  within  the  uncertainty  limits.  The  differences  in  sulfur  conversion  with 
different  sulfur  content  are  most  likely  due  to  systematic  errors  in  the  measurements,  not 
to  a  significant  change  in  sulfur  chemistry  over  the  range  of  sulfur  loadings.  It  can  be 
seen  from  Table  5  that  for  all  three  fuels,  the  SO2  emission  indices  or  sulfur  conversion 
fractions  from  the  three  techniques  agree  within  the  standard  deviations  of  their  data  sets, 
even  without  considering  any  estimates  of  systematic  errors. 

The  overall  result  of  all  three  measurement  techniques  is  that  the  measured  SO2  emission 
represents  the  majority  of  the  fuel  sulfur  content  (see  Table  5).  Measurements  using  the 
TDL  technique  provided  an  upper  limit  for  emitted  S03  concentrations,  as  the  S03 
emissions  were  below  the  detection  sensitivity.  This  upper  limit  is  sufficiently 
constraining  for  the  high  sulfur  fuel  to  verify  the  conclusion  that  most  of  the  emitted 
sulfur  is  in  the  form  of  SO2. 

5.3.  Particulate  Data 

The  UMR  MASS  was  employed  to  characterize  the  particulate  emissions  over  the  range 
of  engine  flight  conditions  and  fuel  types  described  in  this  report.  The  data  characterize 
the  particulates  generated  in  the  engine  and  that  continue  to  exist  12  cm  downstream  of 
the  exit  plane  of  the  engine  exhaust  nozzle.  This  is  the  first  database  of  its  kind  providing 
particulate  characterization  (i.e.  total  concentration,  particulate  number-  and  mass-based 
emission  indices,  soot  volume  fraction,  total  particulate  size  distributions)  as  a  function  of 
engine  operating  parameters  and  fuel  formulation.  A  separate  study  of  JP8+100  fuel  was 
also  performed  and  the  results  of  this  study  are  also  presented  here. A  complete 
description  of  the  measurement  system,  results  and  conclusions  for  the  particulate  data 
are  presented  in  Appendix  D. 


6.  CONCLUSIONS 

A  comprehensive  set  of  data,  including  gaseous,  particulate  precursor,  and  particulate  emissions, 
has  been  acquired  at  simulated  sea-level-static  and  altitude  conditions  for  four  different  fuels. 
Three  JP-8  fuels  with  varying  sulfur  content  from  very  low  (less  than  20  ppmm)  through 
intermediate  (152  ppmm)  to  high  (1113  ppmm)  sulfur  were  tested,  with  the  lowest  sulfur  fuel 
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using  two  sets  of  additives.  An  additional,  advanced  fuel,  JP-8+100  with  a  fuel  sulfur  content  of 
336  ppmm,  was  also  tested.  This  database  provides  valuable  information  for  the  atmospheric 
assessment  of  aerosol  and  aerosol  precursor  emissions  from  aircraft  engines. 

S02  emission  indices  were  measured  by  three  different  techniques,  two  extractive  and  one  non- 
intrusive.  Extractive  samples  were  analyzed  by  ultraviolet  absorption  (gas  analysis  system)  and 
chemical  ionization  mass  spectrometry.  In  situ  measurements  were  performed  with  infrared 
absorption  using  a  tunable  diode  laser.  The  S02  emissions  indices  for  all  three  techniques  agree 
within  their  respective  uncertainty  estimates.  The  measured  S02  emission  represents  the  majority 
of  the  fuel  sulfur  content  for  all  of  the  fuels  used  under  all  engine  and  altitude  conditions. 

Gaseous  emissions  data  were  acquired  by  the  extractive  engine  exhaust  analysis  system  at  all  test 
conditions.  In  part  of  the  testing  period,  the  CO2  analyzer  mal-functioned.  On  those  occasions, 
emission  indices  were  calculated  by  O2  measurements.  The  emission  indices  of  NOx  have 
stronger  dependency  on  combustor  inlet  temperature  (T3)  and  pressure  (P3),  weaker  dependency 
on  fuel  flow  and  no  dependency  on  fuel  sulfur  content.  No  SO2  data  were  obtained  during  the 
low  sulfur  fuel  test  since  the  concentrations  are  below  the  detection  limit  of  this  system. 
Emission  indices  of  SO2  showed  no  dependency  on  T3. 

The  CIMS  system  provided  a  comprehensive  database  on  S02  for  all  fuels.  The  emission  indices 
for  S02  are  essentially  independent  of  engine  exhaust  temperature  and  altitude.  It  is  the  first  time 
that  a  HNO,  database  has  been  obtained  for  engine  emissions  at  the  exit  plane.  Within  the 
uncertainties  in  the  data,  the  emission  indices  for  HNO,  are  independent  of  the  sulfur  level  in  the 
fuel.  The  HNO,  level  decreases  with  temperature,  presumably  because  HNO,  is  not  stable  at  high 
temperatures.  No  CIMS  measurements  of  SO,  were  obtained,  due  to  losses  in  the  sampling  line. 

The  TDL  system  performed  reliably  throughout  the  test  period,  making  measurements  in  both 
double-pass  and  multi-pass  configurations.  During  most  of  the  test,  S02  and  SO,  spectral  regions 
were  measured,  each  of  which  also  contained  several  H,0  lines  which  were  used  as  reference 
lines  for  emission  index  calculations.  For  the  lowest  sulfur  fuel  case,  when  the  S02  levels  were 
below  the  detection  sensitivity,  NO  and  C02  absorption  was  measured.  S02  emission  indices 
were  obtained  for  all  engine  power  conditions  at  altitudes  of  9,144m  (30,000  ft)  and  above  for  all 
but  the  lowest  sulfur  fuel.  SO,  was  always  below  the  detection  sensitivity  of  the  instrument  but, 
for  the  high  sulfur  fuel  case,  the  measurements  allow  an  upper  bound  for  the  SO,  level  to  be 
obtained.  This  upper  bound  is  consistent  with  the  level  of  sulfur  emission  present  as  S02. 

The  MASS  system  provided  a  comprehensive  database  on  particulates.  For  the  total 
concentration,  no  strong  dependence  of  emission  index  on  T3  or  altitude  is  observed.  In  some 
cases  a  weak  dependence  is  observed  where  the  El  appears  to  peak  at  Tj’s  around  600-700  °K. 
The  mean  emission  indices  for  the  high  and  medium  sulfur  cases  are  comparable.  The  mean 
emission  index  for  the  low  sulfur  case  studied  was  significantly  lower  (by  a  factor  of  between  3 
to  4)  than  that  for  the  high  and  medium  sulfur  cases.  However,  each  fuel  had  a  different  additive 
package.  The  mean  number-based  emission  index  for  the  JP8+100  fuel  case  was  (2.0  ±  0.4)  x 
1014  particles/kg  fuel  burned.  This  result  is  comparable  to  the  emission  indices  measured  for  the 
medium  and  high  sulfur  cases.  In  this  case,  however,  yet  a  different  additive  package  is  present 
in  a  different  baseline  fuel. 
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An  analysis  of  the  particle  mass  distributions  indicated  >90%  of  the  particle  mass  fell  within  the 
measurement  regime  of  the  DMA  (10  -  300  nm)  and  thus  the  estimation  of  mass-based  emission 
indices  from  the  mean  volume  diameter  of  the  distributions  and  an  assumed  density  for  the 
particles  was  valid. 

For  all  fuels  and  test  conditions  studied  the  size  distributions  are  of  a  log-normal  type  with  mean 
diameters  in  the  range  50-65  nm.  In  the  case  of  the  medium  sulfur  fuel  study,  data  were  acquired 
for  the  same  engine  operating  conditions  at  different  times  when  test  conditions  were  repeated 
with  the  same  fuel  on  different  days.These  data  provide  an  opportunity  to  explore  the  engine 
performance  issue  of  particulate  emission  stability  between  subsequent  measurements  at 
equivalent  test  conditions.  The  emissions  reproducibility  is  good.  The  standard  deviation  in  the 
mean  diameters  was  1.5%  and  that  of  the  half  widths  of  the  distribution  was  3%. 
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Table  1. — Wright-Patterson  Laboratory  Fuel  Analysis 


TEST 

METHOD 

LIMITS 

MIN 

LIMITS 

MAX 

Sample 

#1 

Sample 

#2 

Sample 

#3 

Sample 

#4 

Sample 

#5 

Total  Acid  Number,  mg  KOH/g 

D974 

0.1 

0.0 

0.0 

0.0 

0.0 

0.008 

Aromatics,  %  vol 

D1319 

22 

15 

17 

10 

15 

9.7 

Aromatics,  % 

HPLC 

14.9 

20.8 

18.2 

20.1 

14.1 

Olefins,  %vol 

5.0 

1.0 

Mercaptan  Sulfur,  %  wt 

D3227 

0.003 

0.000 

0.000 

0.000 

0.000 

0.000 

Total  Sulfur,  %  wt 

D4294 

0.3 

0.0 

0.0 

0.1 

0.0 

0.03 

Sulfur  Content,  ppm 

ASTM3120 

<1 

Sulfur  Content,  ppm 

x-ray 

fluorescence 

~0 

m 

1113 

152 

336 

Distillation 

D86 

IBP,  °C 

Report 

161 

10%  recovered,  °C 

205 

161 

162 

184 

162 

176 

20%  recovered,  °C 

182 

50%  Recovered,  °C 

Report 

179 

180 

214 

184 

197 

90%  recovered,  °C 

Report 

228 

229 

243 

236 

234 

FBP,  °C 

300 

270 

273 

270 

272 

263 

Residue,  %  vol 

1.5 

1.2 

1.2 

1.2 

1.2 

1.2 

Loss,  %  vol 

1.5 

1.3 

1.4 

1.1 

i.i 

1.3 

Flash  Point,  °C 

D56 

38 

39 

38 

51 

40 

45 

Density  at  15  °C,  kg/cu  meter 

D1298 

775 

840 

791 

792 

812 

795 

46 

Freezing  Point,  IC 

D2386 

-40 

-52 

-5l1 

-41 

-51 

-56 

Viscosity  @  -20  °C  cSt 

D445 

8 

3.1 

3.1 

5.7 

3.4 

3.3 

Hydorgen  Content,  %wt 

D3343 

13.4 

14.0 

13.9 

14.0 

14.0 

14.2 

Smoke  Point,  mm 

D1322 

25 

24## 

24## 

24## 

24# 

24.0## 

Copper  Corrosion 

D130 

1 

la 

la 

la 

la 

la 

Thermal  Stability 

D3241 

Tube  Rating  Visual 

<3 

i 

4## 

i 

4## 

i 

Change  in  Press.,  mm  of  Hg 

25 

0 

1 

0 

3 

0 

Existent  Gum,  mg/  IOOmL 

D381 

7 

3 

3 

2 

2 

3 

Water  Reaction 

D1094 

IB 

1 

1 

2# 

1 

2## 

Conductivity,  pS/m 

D2624 

200 

600 

0## 

892## 

0## 

860## 

450 

FSII 

D5006 

0.10 

0.15 

0.00## 

0.11 

0.00## 

0.14 

0.13 

D3338 

42.8 

43.4 

43.3 

43.3 

43.3 

43.4 

Lubricity  (BOCLE) 

D5001 

Report 

0.55 

0.50 

0.65 

0.5 

0.55 

##  — out  of  range 
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Table  2. — Martel  Laboratory  Fuel  Analysis 


Analytical  Parameter 


Filtration  Time 


Particulate  Contamination 


Appearance  /  Workmanship 


Water  Reaction 


Water  Reaction  -  Interface 


Water  Reaction  -  Separation 


Conductivity  of  Aviation  Fuels 


Hydrogen  Content  of  Aviation  Fuels 


Determination  of  Fuel  Icing  Inhibitors 


Hydrocarbon  Constituents  by  F1A: 


Aromatics 


Olefins 


Saturates 


Smoke  Point  of  Aviation  Fuels 


Jet  Fuel  Thermal  Oxidation  Test  (JFTOT) 


Existent  Gum  in  Fuels  by  Jet  Evaporation 


Corrosion  (Copper  Strip  Test) 


Distillation  of  Petroleum  Products 


Initial  Boiling  Point 


10%  Recovered 


20%  Recovered 


30%  Recovered 


40%  Recovered 


50%  Recovered 


60%  Recovered 


70%  Recovered 


80%  Recovered 


90%  Recovered 


End  Point 


Recovery 


Residue 


Loss 


Flash  Point  (Pensky-Martens  closed  cup) 


Freezing  Point 


Gravity,  API  (Hydrometer) _ 


Kinematic  Viscosity  @  -20  deg  C 


Heat  of  Combustion  (Lower  Heating  Value) 


Total  Acid  Number 


Mercaptan  Sulfur 


Sulfur  Non -dispersive  X-Ray  Spectromet 


Sulfur 


Iron 


Copper 


Zinc 


Method 


IASTM  SPEC 


ASTM  D5452 


ASTM  D4176 


ASTM  D1094 


ASTM  D 1 094 


ASTM  D1094 


IASTM  D1 319 


ASTM  D1319 


ASTM  D1319 


ASTM  D1319 


ASTM  D1322 


ASTM  D3241 


ASTM  D3241 


IASTM  D381 


ASTM  D130 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


IASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


ASTM  D86 


|ASTM  D93 


ASTM  D2386 


ASTM  D287 


IASTM  D445 


|ASTM  D240 


ASTM  D664 


ASTM  D3227 


ASTM  D4294 


|ASTM  D4951 


EPA  200.7 


EPA  200.7 


EPA  200.7 


#1 


33.3 


0.2 


bright 

clear 


#2  #3 


33.1  48.7 


0.1  0.4 


yellow 


ASTM  D2624 

106 

850 

mm 

670 

402 

ASTM  D3343 

13.94 

13.89 

13.82 

13.86 

14.02 

ASTM  D5006 

<0.01 

0.11 

0.11 

0.12 

18690 


0.01 


<0.001 


<0.01 


<0.01 


<1 


<1 


<10 
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Table  3. — Fuel  Additives 


Product 

Company 

Quantity 

Components 

Corrosion 

Inhibitor 

Octel 

America 

9  -  22.5 
g/m3  fuel 

70-80%  NJ  Trade  Secret  Registry  #00850201 00 1-5000P, 
20-30%  Xylene,  0-5%  Ethylbenzene,  <250ppm  Benzene 

Conductivity 

Additive 

Octel 

America 

1.5  ppm 
by  weight 

50-60%  Toluene,  5-10%  Mixed  Aromatic  Solvents  C9-16, 
0-5%  Isopropyl  Alcohol,  1-10%  NJ  Trade  Secret  Registry  # 
00850201 00 1-5457P,  1-10%  Dodecyl  Benzene  Sulfonic 
Acid,  10-20%  NJ  Trade  Secret  Registry  #00850201001- 
5037P,  <595ppm  Benzene 

Icing 

Inhibitor 

0.12%  by 
volume 

Diethylene  Glycol  Methyl  Ether 
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Table  4. — Test  Matrix 


Altitude 

t3 

9,144m 

(30,000ft) 

443°C  (830°F) 

SLS 

idle  (~224°C/435°F), 

343°C(650°F),  399°C(750°F), 
443°C(830°F),  mil(~522°C/972°F) 

9,144m 

(30,000ft) 

idle(~226°C/435°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~486°C/907°F) 

12,192m 

(40,000ft) 

idle(~260°C/500°F),  343°C(650°F), 
399°C(750°F),  mil(~447°C/836°F) 

15,240m 

(50,000ft) 

idle(~349°C/660°F),  399°C(750°F), 
mil(~438°C/82 1  °F) 

16,764m 

(55,OOOft) 

idle(~4 1 0°C/769°F), 
mil(~451°C/843°F) 

SLS 

idle  (~218°C/425°F),  ] 

343°C(650°F),  399°C(750°F), 
443°C(830°F),  mil(~5 1 5°C/959°F)  < 

9,144m 

(30,000ft) 

idle(~228°C/443°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~486°C/907°F) 

12,192m 

(40,000ft) 

idle(~264°C/507°F),  343°C(650°F), 
399°C(750°F),  mil(~437°C/818°F) 

15,240m 

(50,000ft) 

idle(~354°C/670°F),  399°C(750°F), 
mi  1(  ~443°C/829°F ) 

16,764m 

(55,000ft) 

mil(~442°C/828°F) 

SLS 

idle  (~221°C/430°F), 

288oC(550°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~5 16°C/960°F) 

9,144m 

(30,000ft) 

idle(~228°C/443°F) ,  343°C(650°F) , 
399°C(750°F),  443°C(830°F), 
mi  1(  ~488°C/9 1 0°F) 

10,668m 

(35,OOOft) 

idle(~228°C/443°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~447°C/837°F) 

12,192m 

(40,000ft) 

idle(~262°C/504°F) ,  343°C(650°F) , 
399°C(750°F),  mil(~442°C/827°F) 

Icing  Inhibitor 


High  sulfur  Jet  A 

( 1 1 13ppm)  with  Corrosion 

Inhibitor 
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Table  4. — Test  Matrix  (Continued) 


15,240m 

(50,000ft) 

idle(~357°C/674°F),  399°C(750°F), 
mil(~443°C/829°F) 

16,764m 

(55,000ft) 

idle(  ~423°C/775°F) , 
mil(~446°C/834°F) 

SLS 

idle  (~227°C/440°F), 

343°C(650°F),  399°C(750°F),  l 

443°C(830°F),  mil(~515°C/959°F)  i 

9,144m 

(30,000ft) 

idle(~228°C/443°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~492°C/9 1 7°F) 

10,668m 

(35,000ft) 

idle(~229°C/445°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~452°C/846°F) 

12,192m 

(40,000ft) 

idle(~269°C/5 17°F),  343oC(650°F), 
399°C(750°F),  mil(~444°C/831°F) 

15,240m 

(50,000ft) 

idle(~355°C/670°F),  399°C(750°F), 
mil(~447°C/838°F) 

16,764m 

(55,000ft) 

idle(~42 1  °C/7 90°F) , 
mil(~45 1  °C/843°F) 

SLS 

idle  (~221°C/430°F), 

343°C(650°F),  399°C(750°F), 
443°C(830°F),  mil(~5 14°C/958°F) 

9,144m 

(30,000ft) 

idle(~233°C/452°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~497°C/926°F) 

10,668m 

(35,000ft) 

idle(~233°C/451°F),  343°C(650°F), 
399°C(750°F),  443°C(830°F), 
mil(~453°C/848°F) 

12,192m 

(40,000ft) 

idle(~269°C/5 17°F),  343oC(650°F), 
399°C(750°F),  mil(~443°C/830°F) 

13,716m 

(45,000ft) 

idle(~3 10°C/590°F),  399°C(750°F), 
mil(~445°C/833°F) 

15,240m 

(50,000ft) 

idle(~352°C/666°F),  399°C(750°F), 
mil(~442°C/828°F) 

16,764m 

(55,OOOft) 

idle(~417°C/782°F), 
mil(  ~454°C/849°F) 

JP8+100  (336ppm) 
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Table  5. — Sulfur  dioxide  Emission  Indices  and  Exhaust  S02  Fractions  of  Fuel  Sulfur, 
Averaged  Over  All  Altitudes  and  Power  Settings 


Fuel 

AFRL  CIMS 

ARI  TDL 

High  Sulfur  Fuel 

2.01  !  0.18 

2.49!  0.62 

1.63  !  0.22 

JP-8  +  100 

0.53  !  0.08 

0.57  !  0.09 

0.56!  0.10 

Mixed  Sulfur  Fuel 

0.31  !  0.065 

0.34!  0.07 

0.26!  0.08 

Sulfur  Dioxide  Emission  Indices  Averaged  Over  Altitude/Power  Settings 
(Error  limits  do  not  include  estimates  of  systematic  errors,  but  are  simply  standard  deviations) 


Fuel 

LeRC  Sampling 

AFRL  CIMS 

ARI  TDL 

High  Sulfur  Fuel 

0.90  !  0.08 

1.19  !  0.28 

0.73  !  0.19 

JP-8  +  100 

0.78  !  0.11 

0.85  !  0.13 

0.83  !  0.22 

Mixed  Sulfur  Fuel 

1.02  !  0.21 

1.11  !  0.22 

0.87  !  0.30 

Fraction  of  Fuel  Sulfur  in  Exhaust  in  the  Form  of  Sulfur  Dioxide 
(ARI  error  limits  include  estimates  of  systematic  errors,  added  in  quadrature  to  standard 
deviations  formed  from  repeated  measurements) 
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Figure  1. — PSL4  Facility  of  NASA  Lewis  Research  Center 
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Fuel/Air  Ratio  Calculated  from  Emissions  Analysis 
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Figure  2. — Horizontal  probe  distributions  of  emissions  fuel/air  ratio  ant  total  temperature 


Vertical  Probe  Position,  cm 

Figure  3. — Vertical  probe  distribution  of  emissions  fuel/air  ratio 
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APPENDIX  A 

Sampling  Probe  Rake  System 

Dennis  Dicki 

Dynacs  Engineering  Company,  Inc. 


A.l.  INTRODUCTION 


The  exhaust  gas  sampling  rake  used  to  draw  off  a  sample  was  a  cruciform  type  rake  which  was 
designed  and  fabricated  by  AEDC  and  used  in  the  previous  1995  engine  emissions  program1. 
The  rake  assembly  consists  of  a  solid  strut  extending  across  the  entire  back  of  the  exhaust  nozzle 
on  the  horizontal  centerline,  with  two  additional  struts  located  on  the  vertical  centerline  which 
are  supported  on  the  rake  support  structure  and  extend  into  the  exhaust  plume  to  a  point  just 
above  or  below  the  center  of  the  horizontal  strut.  This  is  shown  in  Figures  A-l  and  A-2.  The 
horizontal  strut  consists  of  alternating  Mach  number/flow  angularity  probes  and  static 
temperature  probes  between  the  gas  sample  probes  which  are  located  on  5.08  cm  centerlines. 
There  are  a  total  of  7  Mach  number/flow  angularity  probes,  10  static  temperature  probes  and  18 
gas  sampling  probes  on  this  rake.  The  tips  of  the  sampling  probes  were  positioned  21.44  cm 
(8.44  in)  downstream  of  the  exhaust  nozzle  exit  plane.  The  bottom  half  of  the  vertical  rake 
contains  9  standard  gas  sampling  probes.  The  top  half  of  the  vertical  rake  contains  9  open-ended 
0.635  cm  O.D.  glass  lined  tubes. 

The  rakes  were  all  mounted  on  a  support  structure  made  of  15.24  cm  square  structural  tubing. 
The  structure  was  positioned  between  the  engine  exhaust  nozzle  and  the  inlet  to  the  exhaust 
collector,  then  welded  to  the  floor.  Also  mounted  on  this  structure  were  the  mirror  tables  used 
for  the  Aerodyne  analyzer. 

The  sample  probe  tips  on  the  rakes  in  the  bottom  vertical  and  both  horizontal  quadrants  are 
shown  in  Figure  A-3  and  were  designed  to  meet  the  SAE  standards  for  emission  instrumentation. 
The  probes  on  the  top  vertical  rake  consisted  of  a  0.635  cm  glass  lined  tube  which  was  pressed 
into  0.794  cm  tube  and  both  were  chamfered  to  an  angle  of  45°.  See  Figure  A-4  for  details.  The 
Mach  Number/flow  angularity  probes  consisted  of  single,  centered  total  pressure  tap,  surrounded 
by  4  static  type  taps  slightly  downstream  of  the  probe  tip,  on  the  tapered  portion  of  the  tip.  This 
probe  was  designed  to  provide  Mach  Number  and  swirl  characteristics  of  the  exhaust  plume 
(Figure  A-5).  The  static  temperature  was  measured  with  thermocouple  probes  as  shown  in  Figure 
A-6. 

The  sample  flowed  through  smooth  bore,  thin  wall  stainless  steel  tubes  between  the  rakes  and  the 
various  pieces  of  equipment.  The  tubes  were  electrically  heated  and  insulated  either  individually 
or  in  bundles  of  3.  Each  individual  assembly  (hose)  had  a  separate  temperature  controller  in 
order  to  maintain  a  minimum  temperature  of  176.7  °C  (350°F).  The  valve  and  pump  box  were 
also  regulated  to  maintain  the  temperature  at  176.7  °C  (350°F). 
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For  the  first  2  nights  of  the  test  program  the  gas  samples  were  drawn  from  the  exhaust  plume 
using  individual  sample  probes,  which  went  from  the  sample  rake  through  individual  electrically 
heated  stainless  steel  tubes  to  a  valve  switching  box,  through  a  heated  pumping  box  and  then  to 
the  analyzers.  The  valve  switching  box  is  a  piece  of  AEDC  hardware  that  has  the  capability  of 
drawing  from  individual  probes  or  ganged  samples  from  the  rake.  It  is  shown  schematically  in 
Figure  A-7.  Individual  samples  were  drawn  from  the  vertical  rake  the  first  night,  then  the 
horizontal  rake  the  second  night  in  order  to  map  the  size  of  the  exhaust  plume.  The  pump  box 
was  used  to  increase  the  sample  gas  pressure  to  some  value  higher  than  atmospheric  which  is 
required  to  get  it  through  the  bank  of  NASA  gas  analyzers. 

For  the  third  and  subsequent  test  nights,  the  sample  extraction  tubing  on  the  horizontal  rake  was 
reconfigured  so  that  the  sample  from  each  of  the  3  inner  most  sample  probes  from  the  engine 
centerline,  were  ganged  together  and  the  sample  sent  through  electrically  heated  stainless  steel 
tubes  to  the  heated  pump  box  then  to  the  LeRC  analyzers.  See  Figure  A-8  for  schematic.  The 
bottom  vertical  rake  was  dedicated  to  UMR,  which  used  the  3  inner  most  gas  sampling  probes. 
The  top  vertical  rake  was  dedicated  to  AFRL.  The  3  innermost  probe  tubes  were  ganged  together 
at  the  base  of  the  rake  and  connected  to  glass  lined,  electrically  heated  hoses  which  brought  the 
sample  outside  of  the  test  cell  to  the  analyzer. 

REFERENCES 

1.  R.  P.  Howard  et  al  “Experimental  Characterization  of  Gas  Turbine  Emissions  at  Simulated 
Flight  Altitude  Conditions” AEDC-TR-96-3,  September  1996. 
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Figure  A-2.— Top-view  of  sampling  probe  rake  and  engine  nozzle  exit 
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STANDARD  SAMPLING  PROBE 


Figure  A-3. — Standard  Sampling  Probe 


Figure  A-4.  — Philips  Sample  Probe 
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Figure  A-5. — Mach/Flow  Angularity  Probe 


Figure  A-6. — Stagnation  Temperature  Probe 
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APPENDIX  B 

Extractive  Gas  Sampling  System  and  Smoke  Meter 

Kurt  Loos  and  Dawn  Noss 
Dynac  Engineering  Company,  Inc. 


B.l.  INTRODUCTION 

Exhaust  gases  drawn  through  the  sampling  rakes  were  diverted  to  a  gas  bench  containing  an 
array  of  analyzers  dedicated  to  examining  the  concentrations  of  specific  gases  contained  in  the 
sample  stream.  Similarly,  smoke  number  values  were  acquired  using  a  standard  smoke  meter. 

B.2.  INSTRUMENTATION 

Extracted  products  of  combustion  were  analyzed  for  levels  of  carbon  dioxide  (CO2),  carbon 
monoxide  (CO),  oxygen  (O2),  nitric  oxide  (NO),  oxides  of  nitrogen  (NOx),  sulfur  dioxide  (SO2), 
and  unburned  hydrocarbons  (THC)  as  detected  by  a  rack  of  gas  analyzers.  Provided  by 
Rosemount  Analytical,  Inc.  of  La  Habra,  California,  the  Rosemount  System  consisted  of  seven 
individual  instruments  feeding  from  a  single  sample  input  line.  Internal  solenoid  valves  tied  into 
a  control  panel  and  programmable  logic  controller  served  to  distribute  the  flow  through  the 
necessary  sample  conditioning  paths  appropriate  for  each  apparatus.  Hot,  wet,  samples  were 
directed  to  both  the  chemiluminescence  analyzers  for  NO  and  NOx  detection  and  the  THC 
analyzer  that  used  flame  ionization  technology.  For  the  remaining  systems,  the  incoming  sample 
was  sent  first  through  a  refrigeration  bath  of  a  glycol  /  water  mixture  to  lower  the  sample  dew 
point  to  below  4.4  °C  (40  °F).  C02  and  CO  levels  were  determined  through  employing  dedicated 
non-dispersive  infrared  systems  while  O2  was  measured  using  a  paramagnetic  detector  system 

Prior  to  each  test  period,  the  Rosemount  System  was  calibrated  through  the  steady  state  data 
acquisition  system,  ESCORT.  Using  a  direct  connection  which  bypassed  the  sample  line, 
nitrogen  gas  was  passed  through  the  entire  gas  bench  for  a  period  of  six  minutes  to  purge  the 
lines  of  any  residual  gases  and  to  serve  as  a  zero  gas  for  the  analyzers.  This  process  was  repeated 
for  all  active  ranges  associated  with  each  apparatus.  Similarly,  span  gases  of  known  species 
concentration  were  introduced  for  an  up  scale  value.  All  gas  concentrations  were  traceable  to  the 
National  Institute  of  Standards  and  Technology  (NIST)  except  for  those  used  in  calibrating  the 
THC  and  SO2  analyzers.  No  NIST  certified  standard  gases  were  obtainable  for  this  purpose.  In 
lieu  of  standard  information,  the  calibration  gas  supplier's  statement  of  concentration  was  used 
for  set-up  purposes. 

Four  sampling  rakes  (see  Appendix  A)  were  each  connected  to  individual  electrically  heated 
sampling  lines  of  0.635  cm  diameter  stainless  steel  maintained  at  160  °C  (320  °F)  and  plumbed 
to  a  heated  valve  and  pump  enclosure  (HVPE).  The  HVPE  contained  directional  valving  and  an 
array  of  five  pumps  (  see  Figure  A-9). 
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To  ensure  that  proper  sample  conditions  were  maintained  for  analyzer  accuracy,  pressure, 
temperature,  and  flow  information  were  monitored  through  ESCORT.  Temperature 
instrumentation  existed  at  the  exit  of  each  sampling  probe,  at  each  junction  of  two  heated  lines, 
throughout  the  HVPE,  and  at  the  input  location  for  both  the  Rosemount  system  and  the  smoke 
meter.  The  amount  of  sample  gas  received  by  the  analyzers  was  displayed  on  meters  within  the 
control  room.  If  the  flow  was  insufficient  for  the  entire  gas  bench  to  be  run  as  a  unit,  sampling 
through  the  Rosemount  System  was  taken  by  an  appropriate  combination  of  individual  analyzers. 


B.3.  DATA  ACQUISITION 

Initially  each  probe  was  sampled  to  determine  whether  the  engine  emissions  profile  would  prove 
to  be  symmetrical.  Once  symmetry  was  determined,  the  four  arms  were  individually  assigned  to 
each  research  group.  This  allowed  parallel  sampling  of  the  emissions.  At  no  time  was  the 
sample  gas  shared  by  any  of  the  analysis  centers. 

Species  level  output  was  sent  from  the  individual  Rosemount  analyzers  to  the  data  system, 
ESCORT,  where  it  was  converted  from  milli-volts  to  appropriate  engineering  units  based  upon 
the  calibration  information  specific  to  that  test  period.  Real  time  updates  of  concentration  values 
were  displayed  once  per  second  for  preliminary  review.  Additionally,  a  data  reduction  routine 
was  incorporated  into  the  ESCORT  programming.  Emission  indices  were  also  calculated 

The  gas  sampling  system  used  was  designed  and  built  by  Rosemount  Analytical,  Inc.  to 
continuously  monitor  aircraft  turbine  engine  exhaust.  The  analyzers  used  in  the  system  and  their 
range  capabilities  were: 


Analyzer 

Operable  Ranges 

Calibration 

Concentrations 

Calibration  Gas 
Constituents 

Analyzer  Technology 

C02-Rosemount  Model  880 

0-2500ppm 

2000ppm 

EBilSIilBi 

Non-Dispersive  IR 

0-5% 

4% 

0-15% 

12.90% 

CO-Rosemount  Model  880 

0-1% 

0.79% 

CO  in  N2 

Non-Dispersive  IR 

0-5% 

2.03% 

0-10% 

8.50% 

02-Rosemount  Model  755 

0-25% 

20.90% 

Air 

Paramagnetic 

NOx-Rosemount  Model  955 

0-100ppm 

92.8ppm 

Chemiluminescence 

O-lOOOppm 

764ppm 

NO-Rosemount  Model  955 

0-100ppm 

92.2ppm 

NO  in  N2 

Chemiluminescence 

;  O-lOOOppm 

763ppm 

SCVRosemount  Model  ETL  9100 

0-10ppm 

7.28ppm 

Electro-Optic  Etalon 

0-100ppm 

THC-Rosemount  Model  404 

0-250ppm 

70ppm 

C3H8  in  Air 

Flame  Ionization 
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The  gas  samples  were  passed  through  heated,  smooth  bore,  stainless  steel  hoses  to  the  sample 
distribution  system.  The  sample  distribution  system  was  previously  used  at  Arnold  Engineering 
Development  Center  for  the  Atmospheric  Effects  of  Aviation  Project,  and  allows  any 
combination  of  probe  samples  to  be  directed  to  the  gas  analyzer  system.  NASA’s  distribution 
system  included  five  individually  controlled  vacuum  pumps  which  raised  the  sample  pressure 
above  atmospheric.  The  pumps  allowed  the  samples  to  be  drawn  into  the  analyzers  from  the  test 
section. 

The  samples  directed  to  the  oxygen,  oxides  of  carbon,  and  sulfur  dioxide  analyzers  were  sent 
through  a  glycol/water  refrigeration  bath  where  condensibles  were  removed  and  the  sample  dew 
point  was  lowered  below  4.44  °C  (40°F).  The  temperature  and  humidity  of  the  samples  directed 
to  the  oxides  of  nitrogen  and  hydrocarbon  analyzers  were  maintained. 

An  SAE  smoke  meter  was  used  as  a  qualitative  measure  of  particles  in  the  exhaust.  The 
measurement  was  reported  as  smoke  number. 
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APPENDIX  C 

NASA  Engine  Parameters  and  Emission  Data 

Changlie  Wey 

Dynacs  Engineering  Company,  Inc. 

Chowen  Chou  Wey 
Army  Research  Lab./NASA  LeRC 


Cl.  INTRODUCTION 

Emissions  measurements  were  performed  in  the  mixed  stream  of  engine  core  flow  and  fan  by¬ 
pass  air  using  a  cruciform  type  sampling  rake  described  in  Appendix  A.  The  tips  of  the  sampling 
probes  were  0.214  meter  behind  the  engine  exhaust  nozzle.  There  are  total  of  36  gas  sampling 
probes,  9  for  each  arm.  Figure  C-l  shows  the  position  of  the  sampling  probes  relative  to  the 
engine  exhaust  nozzle.  Heated  stainless  steel  sampling  tubes  connect  sampling  probes  to  the 
heated  pump  box,  which  is  located  outside  the  chamber.  The  temperature  of  all  the  components 
of  the  sampling  system  downstream  of  the  sampling  probe  was  maintained  at  or  above  176  °C. 

Gaseous  emissions  in  the  engine  exhaust  were  measured  by  an  aircraft  turbine  engine  exhaust 
analysis  system.  The  system,  described  in  Appendix  B,  is  similar  to  a  system  that  would  be  used 
for  commercial  aircraft  engine  emissions  certification  testing  and  therefore  should  provide 
reliable  data  for  this  study.  Operation  of  the  system  followed  the  ICAO  regulations1. 

Concentrations  of  CO2,  CO,  O2,  NO,  NOx,  SO2,  and  total  unburned  hydrocarbon  were  measured. 
Fuel/air  ratio,  combustion  efficiency  and  emission  indices  were  calculated  from  the  measured 
concentrations  per  the  SAE  ARP15332. 

A  smoke  meter  (ROSECO  473B)  was  used  to  measure  smoke  per  ICAO  standard.  The 
reflectances  of  sampled  filter  paper  were  measured  by  a  reflection  densitometer  and  the  result 
was  recorded  as  smoke  number. 

Emissions  data  were  obtained  for  all  test  conditions  with  engine  power  level  ranging  from  idle  to 
takeoff  and  altitudes  ranging  from  sea  level  to  15  Km.  The  spatial  profiles  of  the  emissions  were 
determined  by  drawing  the  gas  samples  from  individual  sample  probes  on  the  first  two  days. 
Based  on  the  distributions,  the  three  innermost  sample  probes  from  the  engine  centerline  of  each 
horizontal  arm  were  ganged  together  for  the  gas  analysis  system  for  the  rest  of  test. 


1  ICAO  International  and  Recommended  Practices  ANNEX  16  as  “Environmental  Protection”,  Volume  II  - 
“Aircraft  Engine  Emissions”,  1993 

2  SAE  Aerospace  Recommended  Practice  ARP1533,  “  Procedure  for  the  Calculation  of  Gaseous  Emissions  from 
Aircraft  Turbine  Engines”,  1996 
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C.2.  RESULT 


All  emission  data  along  with  the  engine  parameters  are  presented  in  Tables  C-l  through  C-5. 
Operation  parameters  listed  in  these  tables  include  simulated  altitude  (ALT),  engine  inlet 
pressure  (PS1),  dew  point  of  inlet  air  (Dewp),  combustor  inlet  pressure  (P3)  and  temperature 
(T3),  engine  air  flow  rate  (WA1),  combustor  air  flow  rate  (WABURN),  fuel  flow  rate  (WFT). 
Fuel/ Air  ratio  is  determined  from  measured  air  flow  and  fuel  flow.  Concentrations  of  gaseous 
species  and  emissions  indices  (El)  are  also  listed. 

NO  and  NOx  emission  indices  are  plotted  in  Figures  C-2  through  C-l  for  respective  altitudes  and 
displayed  as  functions  of  combustor  inlet  temperature,  T3.  Figures  C-8  through  C-l 3  contain  the 
same  information  as  functions  of  combustor  inlet  pressure,  P3.  It  is  obvious  that  NOx  emission 
indices  have  no  dependency  on  fuel  sulfur  contents. 

Figures  C-14  through  C-18  show  the  NOx  emission  indices  as  a  function  of  fuel  flow  for  each 
fuel.  Figures  C-19  and  C-20  show  the  NOx  emission  indices  as  a  function  of  T3  and  P3 
respectively  in  case  of  the  high  sulfur  fuel.  Compare  these  2  plots  with  Figure  C-l 5,  one  can 
conclude  that  NOx  emission  indices  have  stronger  dependency  on  T3  and  P3  than  that  of  fuel 
flow.  Figure  C-21  shows  a  typical  comparison  of  the  calculated  fuel-air  ratio  by  using  emissions 
measurements  and  by  using  fuel  flow  measurements. 

There  are  no  SO2  data  obtained  during  the  low  sulfur  fuel  tests  since  the  concentrations  are 
below  the  detection  limit.  El  (SO2)  values  are  plotted  in  Figures  C-22  through  C-24  for 
respective  fuels  and  displayed  as  functions  of  combustor  inlet  temperature,  T3.  It  shows  no 
dependency  on  T3. 

Figure  C-25  shows  the  smoke  number  versus  altitude  at  a  constant  combustor  inlet  temperature. 
It  shows  that  the  smoke  number  decreases  as  the  altitude  increases.  It  can  also  be  considered  as 
the  smoke  number  increases  as  the  combustor  inlet  pressure  increases. 
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Table  C-l. — Summary  of  Fuel  lf  Test  Result 


Fuel  1:  low  sulfer  fuel  with  corrosion  inhibitor 


Table  C-l. — Summary  of  Fuel  lf  Test  Result  (Continued) 


Fuel  1:  low  sulfer  fuel  with  corrosion  inhibitor 


Table  C-2. — Summary  of  Fuel  2*  Test  Result 


Fuel  2 :  low  sulfur  fuel  with  corrosion  inhibitor,  conductivity  additive  and  icing  inhibitor 


Table  C-2. — Summary  of  Fuel  2*  Test  Result  (Continued) 
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Fuel  2  :  low  sulfur  fuel  with  corrosion  inhibitor,  conductivity  additive  and  icing  inhibitor 


Table  C-3. — Summary  of  Fuel  3*  Test  Result 
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Fuel  4:  mixed  medium  sulfur  fuel  with  corrosion  inhibitor,  conductivity  additive  and  icing  inhibitor 
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Fuel  4:  mixed  medium  sulfur  fuel  with  corrosion  inhibitor,  conductivity  additive  and  icing  inhibitor 


Table  C-5 — Summary  of  Fuel  5*  Test  Result 


Fuel  5:  JP8  + 100 


Table  C-5 — Summary  of  Fuel  5*  Test  Result  (Continued) 
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Table  C-5 — Summary  of  Fuel  5*  Test  Result  (Continued) 
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Fuel  5:  JP8  + 100 
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Figure  C-l. — Sampling  Probe  Locations  relative  to  the  Engine  Exit  Plane 
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Figure  C-2.— EI(NO)  and  EI(NOx)  vs.  T3  at  Sea  Level 


Figure  C-3.— EI(NO)  and  EI(NOx)  vs.  T3  at  9.14  Km 


Figure  C-5.— EI(NO)  and  EI(NO*)  vs.  T3  at  12.19  Km 


24  Km 


Figure  C-7. — EI(NO)  and  EI(NO*)  vs.  T3  at  16.76  Km 
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Figure  C-8. — EI(NO)  and  EI(NOx)  vs.  P3  at  Sea  Level 
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Figure  C-9. — EI(NO)  and  EI(NOx)  vs.  P3  at  9.14  Km 
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Figure  C-10.- 
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Figure  C-ll.— EI(NO)  and  EI(NO*)  vs.  P3  at  12.19  Km 
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>)  and  EI(NOx)  vs.  P3  at  15.24  Km 
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Figure  C-13.— EI(NO)  and  EI(NOx)  vs.  P3  at  16.76  Km 
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Figure  C-14. — EI(NOx)  vs.  fuel  flow  for  Fuel  1 
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Figure  C-15. — EI(NOx)  vs.  fuel  flow  for  Fuel  2 
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Figure  C-16. — EI(NOx)  vs.  fuel  flow  for  Fuel  3 
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Figure  C-17. — EI(NOx)  vs.  fuel  flow  for  Fuel  4 
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•EI(NOx)  vs.  fuel  flow  for  Fuel  5 
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Figure  C-19. — EI(NOx)  vs.  combustor  inlet  temperature  for  Fuel  3 
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combustor  inlet  pressure  for  Fuel  3 
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Figure  C-21. — Comparisons  of  Fuel-Air  Ratio 
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Figure  C-22. — EI(SC>2)  vs.  T3  for  Fuel  3  (high  sulfur  fuel  ~1113ppm) 
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Figure  C-23. — EI(S02>  vs.  T3  for  Fuel  4  (medium  sulfur  fuel  ~152ppm  ) 
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for  Fuel  5  (JP8+100  fuel  ~336ppm  S) 


Figure  C-25, 


APPENDIX  D 

Particulate  Sampling  Measurements 


D.E.  Hagen,  P.D.  Whitefield,  M.B.  Trueblood,  M.E.  Wilson  and  D.  Olson 
Cloud  and  Aerosol  Sciences  Laboratory 
University  of  Missouri-Rolla,  Rolla  MO  65401,  USA 


D.l.  INTRODUCTION 

The  industry  standard  for  particulate  measurements  is  the  smoke  number.  Smoke  number 
although  of  considerable  use  for  engine  design  considerations,  does  not  provide  the 
fundamental  physical  characterization  of  engine  particulate  emissions  required  to  support 
models  that  explore  the  environmental  impact  of  these  emissions  (Paladino  1997).  These 
models  require  inputs  on  such  particulate  physical  characteristics  as:  concentration, 
particulate  number-based  and  mass-based  particulate  emission  index,  size  distribution, 
hydration/growth  potential  (as  described  by  measurements  of  particulate  size  dependent 
soluble  mass  fraction)  and  morphology.  The  University  of  Missouri-Rolla  Mobile 
Aerosol  Sampling  System  (UMRMASS)  has  been  specifically  designed  to  measure  these 
desired  physical  characteristics  and  has  been  employed  successfully  in  previous  altitude 
chamber,  ground  test,  combustor  rig  and  in  flight  test  environments  (see.  for  example 
Howard  et  al  1996,  Hagen  et  al  1992,  1993a,  1996,  1998  and  Paladino  1998).  In  this 
appendix  we  describe  the  application  of  UMRMASS  to  characterize  the  particulate 
emissions  of  an  F100-E200  series  engine  as  a  function  of  combustor  inlet  temperature 
(T3),  altitude  and  fuel  formulation.  This  study  also  included  the  examination  of  the 
particulate  emissions  from  the  same  engine  burning  JP8  with  the  +100  fuel  additive. 


D.2.  EXPERIMENTAL  APPROACH 

A  trailer-based  UMRMASS  was  employed  in  this  NASA  Lewis  sponsored  F100  altitude 
chamber  emissions  study.  The  trailer  was  located  adjacent  to  the  NASA  LeRC  PSL 
facility  and  at  a  distance  of  approximately  30m  from  the  engine  test  cell.  Engine  exhaust, 
extracted  with  the  sampling  probe  rake  was  diluted  with  heated  dry  particle  free  air 
immediately  downstream  of  the  probe  manifold  assembly  and  the  resulting  dried  and 
diluted  exhaust  sample  was  ducted  to  the  trailer  through  a  heated  stainless  steel  sample 
line  (0.0064  m  o.d.  at  423  °K).  A  schematic  diagram  of  the  UMRMASS  configuration  in 
the  trailer  is  given  in  Figure  D-l.  The  diluted  sample  arriving  at  the  trailer  was  distributed 
to  the  following  measurement  stations: 

(1)  Total  concentration  -  the  total  concentration  was  recorded  in  real  time  using  a 
commercially  available  condensation  nucleus  counter  (CNC)  with  a  1Hz  sample 
frequency.  Since  CNC’s  will  saturate  at  high  particulate  concentrations  (Hagen  et  al. 
1993b)  farther  particulate  dilution  in  the  total  concentration  measurement  line  was 
achieved  using  needle  filter  diluters  (Olson  et  al.  1996).  The  total  concentration  was 
monitored  continuously  throughout  the  test  and  was  recorded  as  a  function  of  time. 
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(2)  Size  distribution  -  size  distributions  were  recorded  using  the  differential  mobility 
analysis  technique  (Howard  et  al.  1996).  Diluted  exhaust  sample  was  either  ducted 
directly  to  the  differential  mobility  analyzer  (DMA)  and  a  size  spectrum  was  recorded 
in  real-time  or  the  diluted  sample  was  ducted  to  a  40  liter  storage  tank  after  which 
size  spectra  were  recorded  off-line  but  typically  within  a  few  minutes  of  storage. 

Since  the  MASS  methodology  has  been  described  extensively  elsewhere  (Hagen  1993a, 
Howard  et  al.  1996)  no  further  discussion  will  be  provided  here. 


D.3.  PARTICULATE  CHARACTERIZATION  RESULTS 

In  this  project  the  particulate  physical  characteristics,  described  in  section  D-2  were 
measured  as  a  function  of  altitude  and  combustor  inlet  temperature  (T3)  using  Jet  A  jet 
fuel  with  3  sulfur  contents  defined  as  low,  medium  and  high  each  including  several 
different  additives,  and  with  JP8  with  the  +100  additive  (see  Tables  1-4).  The  data  are  all 
reported  using  the  escort  number  for  reference  to  the  test  matrix.  The  data  presented 
below  is  based  on  the  entire  test  program  data  set  excluding  test  point  periods  406-436 
and  450-498. 

Tables  D-l  thru  D-4  are  compilations  of  the  mean  particulate  number-based  El’s,  mean 
mass-based  El’s,  and  particulate  volume  fractions  for  each  test  point  with  each  table 
representing  the  data  for  a  given  fuel.  The  quoted  uncertainty  in  all  cases  is  based  on  the 
analysis  of  both  statistical  and  systematic  errors  arising  from  the  measurement 
methodologies.  Table  D-5  lists  the  mean  particulate  number-based  El’s,  mean  mass- 
based  El’s,  and  particulate  volume  fractions  for  all  measurements  on  a  given  fuel.  The 
mean  mass-based  El’s  and  particulate  volume  fractions  are  calculated  using  the  measured 
size  distribution  and  number-based  El’s  for  each  test  condition.  An  analysis  of  the 
particle  mass  distributions  indicated  >90%  of  the  particle  mass  fell  within  the 
measurement  regime  of  the  DMA  (10-300nm),  see  Figures  D-2  and  D-3.  Figures  D-4 
through  D-7  plot  the  number-based  El’s  for  all  fuels  as  a  function  of  T3  for  fixed 
altitudes.  Associated  with  each  escort  number  there  are  mean  size  distributions.  These 
size  distributions  are  plotted  in  Figures  D-8  thru  D-35.  Table  D-6  reports  the  average 
mean  diameters  for  all  size  distributions  measured  at  all  engine  operating  conditions  for 
each  fuel  studied.  The  quoted  uncertainty  is  one  standard  deviation.  Where  more  than  one 
size  distributions  can  be  associated  with  a  given  test  condition  (as  a  result  of  returning  at 
a  later  time  to  any  given  test  condition)  each  distribution  associated  with  the  same  test 
condition,  and  identified  by  its  escort  number  is  plotted  on  the  same  axis  in  Figure  D-36. 
These  data  are  presented  to  explore  issues  of  the  engines  reproducibility  of  particulate 
emission  characteristics  for  a  given  set  of  operational  parameters. 
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D.4.  CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  particulate  emissions  measurements 

presented  in  section  D-3. 

Emission  Indices: 

•  No  strong  dependence  of  El  on  T3  or  altitude  is  observed.  In  some  cases  a  weak 
dependence  is  observed  where  the  El  appears  to  peak  at  Ts’s  around  600-700  °K  (see 
Figures  D-4  -  D-7). 

•  The  lack  of  any  strong  dependence  of  El  on  T3  or  altitude  justifies  the  representation 
of  the  fuel  effect  on  particulate  emissions  using  the  mean  value  for  all  test  conditions 
studied  for  a  given  fuel  (Table  D-5). 

•  The  mean  Els  for  the  high  and  medium  sulfur  cases  are  comparable.  The  mean  El  for 
the  low  sulfur  case  studied  was  significantly  lower  (by  a  factor  of  between  3  to  4) 
than  that  for  the  high  and  medium  sulfur  cases.  There  were,  however,  several 
variables  in  the  fuel  formulation  beyond  that  of  sulfur  content.  The  base  fuel  for  the 
low  sulfur  case  was  hydrotreated  and  contained  a  different  additive  package  to  the 
high  sulfur  which  had  not  been  hydrotreated.  The  medium  sulfur  fuel  was  a  blend  of 
the  high  and  low  sulfur  fuels  providing  the  desired  mid-range  sulfur  concentration.  It 
is  clear  from  these  data  that  particulate  emissions  can  be  significantly  affected 
through  fuel  formulation  modification  but  it  is  not  clear  what  facet  of  modification  is 
responsible  for  observed  emissions  changes. 

•  The  mean  number  -based  El  for  the  JP8+100  fuel  case  was  (2.0  ±  0.4)  x  1014 
particles/kg  fuel  burned.  This  result  is  comparable  to  the  Els  measured  for  the 
medium  and  high  sulfur  cases.  In  this  case,  however,  yet  a  different  additive  package 
is  present  in  a  different  baseline  fuel. 

•  The  mean  mass-based  El’s  and  particulate  volume  fractions,  both  of  which  are 
calculated  using  the  measured  size  distribution  and  number-based  El’s  for  each  test 
condition,  exhibit  the  same  trends  as  a  function  of  fuel  formulation. 

•  An  analysis  of  the  particle  mass  distributions  indicated  >90%  of  the  particle  mass  fell 
within  the  measurement  regime  of  the  DMA  (10  -  300  nm)  and  thus  the  estimation  of 
mass-based  El’s  from  the  mean  volume  diameter  of  the  distributions  and  an  assumed 
density  for  the  particles  was  valid. 

Size  Distributions 

•  For  all  fuels  and  test  conditions  studied  the  size  distributions  are  of  a  log-normal  type 
with  mean  diameters  in  the  range  50-65nm  (see  Table  D-6). 

In  the  case  of  the  medium  sulfur  fuel  study,  data  were  acquired  for  the  same  engine 

operating  conditions  at  different  times  when  test  conditions  were  repeated  with  the  same 
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fuel  on  different  days.  These  data  provide  an  opportunity  to  explore  the  engine 
performance  issue  of  particulate  emission  stability  between  subsequent  measurements  at 
equivalent  test  conditions.  The  emissions  reproducibility  is  good.  The  standard  deviation 
in  the  mean  diameters  was  1.5%  and  that  of  the  half  widths  of  the  distribution  was  3%. 
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ESP* 

PTE 

mm 

Fuel 

El  (n)  * 

El(n)  unc. 

buss 

El(m)  unc. 

PVF* 

PVF  unc. 

596 

El 

511 

MSE IB 

1.03E+14 

+/-  2.07E+13 

0.01 

+/-  2.77E-03 

1.66E-11 

+/-  3.32E-12 

599 

O 

616 

WEEM 

2.99E+14 

+/-  5.99E+13 

0.11 

+/-  2.22E-02 

5.07E-10 

+/-  1.01  E-10 

603 

B 

672 

nm 

2.62E+14 

+/-  5.25E+13 

0.10 

+/-  1.94E-02 

6.31  E-10 

+/-  1.26E-10 

607 

B 

716 

EH 

1.02E+14 

+/-  2.04E+13 

0.08 

+/-  1.60E-02 

6.98E-10 

+/-  1.40E-10 

611 

B 

755 

Ena 

1.61E+14 

+/-  3.22E+13 

0.25 

+/-  5.00E-02 

2.73E-09 

+/-  5.46E-10 

569 

m 

537 

WENM 

1.59E+14 

+/-  3.18E+13 

0.03 

+/-  5.43E-03 

3.20E-11 

+/-  6.41  E-1 2 

572 

m 

616 

WEBtSM 

4.79E+14 

+/-  9.59E+13 

0.13 

+/-  2.69E-02 

3.09E-10 

+/-  6.18E-11 

575 

EEl 

671 

ksb 

5.66E+14 

+/-  1.13E+14 

0.32 

+/-  6.44E-02 

9.99E-10 

+/-  2.00E-10 

580 

EE1 

716 

WEFM 

4.69E+14 

+/-  9.38E+13 

0.26 

+/-  5.22E-02 

1.1  IE-09 

+/-  2.23E-10 

584 

EH 

628 

■SMS 

1.25E+14 

+/-  2.50E+13 

0.05 

+/-  9.83E-03 

6.48E-1 1 

+/-  1.30E-11 

587 

EH 

672 

EH9 

2.04E+14 

+/-  4.08E+13 

0.06 

+/-  1.12E-02 

9.21  E-11 

+/-  1.84E-11 

591 

EH 

716 

U1MM 

2.12E+14 

+/-  4.25E+13 

0.11 

+/-  2.10E-02 

2.32E-10 

+/-  4.64E-1 1 

594 

m 

716 

mssm 

1.24E+14 

+/-  2.47E+13 

0.06 

+/-  1.21E-02 

9.72E-1 1 

+/-  1.94E-11 

ESP  Escort  parameter  /  number 
Alt  Altitude 

El  (n)  Number-based  emission  index 
El  (m)  Mass-based  emission  index 
PVF  Particle  volume  fraction 
unc  uncertainty  (Is). 

Table  D-l:  Emission  Indices  and  Soot  Volume  Fractions  as  a  Function  of 
Altitude,  T3  and  Fuel  -  High  Sulfur 
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ESP* 

Alt* 

T3 

Fuel 

El  (n)  * 

El(n)  unc. 

El(m)* 

El(m)  unc. 

PVF* 

PVF  unc. 
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+/-  2.98E-10 
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El 
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+/-  3.53E+13 

0.14 

1.16E-09 

+/-  2.33E-10 

835 

O 

616 

JP8+100 

2.85E+14 

+/-  5.69E+13 

0.10 

+/-  2.02E-02 

2.15E-10 

+/-  4.29E-11 

837 

O 

672 

JP8+100 

1.76E+14 

+/-  3.53E+13 

0.08 

+/-  1.65E-02 

4.24E-10 

+/-  8.49E-11 

840 

El 

Mil 

JP8+100 

1.3E+14 

+/-  2.60E+13 

0.07 

+/-  1.43E-02 

5.07E-10 

+/-  1.01E-10 

844 

O 

766 

JP8+100 

6.02E+13 

+/-  1.20E+13 

0.03 

+/-  6.67E-03 

3.63E-10 

+/-  7.26E-11 

890 

in 

505 

JP8+100 

3.48E+13 

+/-  6.97E+12 

0.00 

+/-  6.47E-04 

6.36E-12 

+/-  1.27E-12 

894 

m 

JP8+100 

4.14E+14 

+/-  8.28E+13 

0.17 

+/-  3.37E-02 

3.97E-10 

+/-  7.94E-11 

899 

ED 

728 

JP8+100 

1.39E+14 

+/-  2.77E+13 

0.07 

+/-  1.42E-02 

5.03E-10 

+/-  1.01E-10 

803 

m 

616 

JP8+100 

4.81  E+14 

+/-  9.63E+13 

0.16 

+/-  3.16E-02 

2.16E-10 

+/-  4.33E-11 

806 

m 

m.rt'M 

JP8+100 

3.09E+14 

+/-  6.19E+13 

0.16 

+/-  3.22E-02 

3.98E-10 

+/-  7.96E-11 

811 

m 

716 

JP8+100 

2.4E+14 

+/-  4.81E+13 

0.14 

+/-  2.83E-02 

5.50E-10 

+/-  1.10E-10 

871 

Eia 

JP8+100 

6.36E+13 

+/-  1.27E+13 

0.01 

+/-  2.34E-03 

3.00E-11 

+/-  6.00E-12 

875 

EH 

678 

JP8+100 

3.19E+14 

+/-  6.38E+13 

0.11 

+/-  2.10E-02 

1.55E-10 

+/-  3.10E-11 

879 

EH 

718 

JP8+100 

2.59E+14 

+/-  5.19E+13 

0.12 

+/-  2.41  E-02 

2.63E-10 

+/-  5.26E-11 

814 

EH 

WiVrM 

JP8+100 

1.4E+14 

+/-  2.80E+13 

0.04 

+/-  8.79E-03 

5.20E-11 

+/-  1.04E-11 

817 

EH 

672 

JP8+100 

2.26E+14 

+/-  4.53E+13 

0.09 

+/-  1.82E-02 

1.31E-10 

+/-  2.63E-11 

822 

EH 

716 

JP8+100 

1.92E+14 

+/-  3.84E+13 

0.10 

+/-  1.94E-02 

1.97E-10 

+/-  3.94E-11 

824 

EH 

689 

JP8+100 

6.92E+13 

+/-  1.38E+13 

0.02 

+/-  4.66E-03 

3.65E-11 

+/-  7.30E-12 

827 

EH 

728 

H235ES 

1.38E+14 

+/-  2.77E+13 

0.07 

+/-  1.36E-02 

8.50E-11 

+/-  1.70E-11 

882 

EH 

689 

JP8+10d 

5.53E+13 

+/-  1.11E+13 

0.02 

+/-  3.31  E-03 

2.74E-1 1 

+/-  5.48E-12 

886 

EH 

722 

JP8+100 

1.18E+14 

+/-  2.36E+13 

0.05 

+/-  1.05E-02 

8.05E-1 1 

+/-  1.61E-11 

* 

ESP  Escort  parameter  /  number 
Alt  Altitude 

El  (n)  Number-based  emission  index 
El  (m)  Mass-based  emission  index 
PVF  Particle  volume  fraction 
unc  uncertainty  (Is) 

Table  D-2:  Emission  Indices  and  Soot  Volume  Fractions  as  a  Function  of 
Altitude,  T3  and  Fuel  -  JP8+100 
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ESP* 

Alt* 

Ts 

Fuel 

El  (n)  * 

El(n)  unc. 

El(m)* 

El(m)  unc. 

PVF  * 

PVF  unc. 

621 

O 

494 

Mids" 

1.80E+13 

+/-  3.59E+12 

0.01 

+/-  2.81  E-03 

5.08E-11 

+/-  1.02E-11 

624 

El 

eh 

1H 

5.76E+13 

+/-  1.15E+13 

0.01 

+/-  2.93E-03 

8.40E-11 

627 

El 

616 

Mid  S 

4.10E+14 

4-/-8.21E+13 

0.24 

+/-  4.70E-02 

1.01E-09 

+/-  2.02E-10 

630 

El 

672 

MidS 

3.86E+14 

MEE33SH 

0.26 

1.37E-09 

+/-  2.74E-10 

633 

El 

716 

MidS 

2.74E+14 

+■/-  5.47E+13 

0.22 

+/-  4.36E-02 

1.76E-09 

+/-  3.52E-10 

636 

El 

789 

MidS 

1.57E+14 

+7-  3.13E+13 

0.13 

+/-  2.62E-02 

1.29E-09 

+/-  2.58E-10 

WH7M 

El 

505 

MidS 

6.56E+13 

BEEPSEI 

E«MI 

+/-  1.26E-03 

1.10E-11 

+/-  2.21  E-12 

711 

O 

616 

MidS 

1.42E+14 

+/-  2.83E+13 

0.07 

+/-  1.47E-02 

3.88E-10 

+/-  7.75E-11 

712 

El 

672 

MidS 

1.93E+14 

+/-  3.87E+13 

0.09 

+/-  1.72E-02 

5.10E-10 

+/-  1.02E-10 

716 

El 

716 

MidS 

9.84E+13 

+/-  1.97E+13 

0.06 

+/-  1.12E-02 

4.99E-10 

+/-  9.98E-11 

718 

El 

761 

MidS 

7.94E+13 

+•/•  1.59E+13 

0.05 

5.50E-10 

+/-  1.10E-10 

720 

El 

561 

MidS 

2.71  E+1 3 

+/-  5.42E+12 

0.01 

+/-  1.44E-03 

4.38E-11 

+/-  8.77E-12 

773 

m 

MidS 

1.60E+14 

BEEE3EI 

HHM 

+/-  1.03E-02 

6.11E-11 

+/-  1.22E-11 

650 

m 

622 

MidS 

1.58E+14 

+/-  3.15E+13 

0.07 

+/-  1.49E-02 

9.71  E-11 

+/-  1.94E-11 

687 

m 

628 

MidS 

1.97E+14 

+/-  3.93E+13 

0.07 

+/-  1.35E-02 

1.18E-10 

+/-  2.36E-11 

653 

EH 

672 

MidS 

1.58E+14 

+/-  3.15E+13 

0.07 

+/-  1.49E-02 

wssssa 

+/-  1.94E-11 

656 

m 

716 

MidS 

3.02E+14 

+•/-  6.04E+13 

0.19 

+/-  3.82E-02 

4.01  E-10 

+/-  8.02E-11 

659 

eh 

686 

MidS 

1.80E+14 

+/-  3.59E+13 

0.09 

+/-  1.72E-02 

9.05E-11 

+/- 1.81  E-11 

784 

m 

700 

MidS 

7.20E+13 

+/-  1.44E+13 

E«I»U 

+/-  5.75E-03 

4.69E-11 

+/-  9.38E-12 

790 

m 

728 

MidS 

1.23E+14 

+/-  2.46E+13 

0.07 

+/-  1.32E-02 

9.76E-11 

+/-  1.95E-11 

640 

BO 

505 

MidS 

9.18E+13 

+/-  1.84E+13 

0.02 

+/-  4.15E-03 

EB3BI 

+/-  5.47E-12 

m 

672 

MidS 

3.23E+14 

+-/-  6.46E+13 

0.17 

+/-  3.33E-02 

KS^d=BIsl 

+/-  8.52E-11 

700 

m 

683 

MidS 

1.88E+14 

+•/•  3.76E+13 

0.06 

+/-  1.26E-02 

6.39E-11 

+/-  1.28E-11 

768 

m 

728 

MidS 

2.39E+14 

+/-  4.78E+13 

0.13 

!KSFi13S^ 

4.74E-10 

+/-  9.49E-11 

646 

eh 

755 

MidS 

2.15E+14 

+/-  4.29E+13 

0.13 

+/-  2.62E-02 

8.69E-10 

+/-  1.74E-10 

695 

eh 

716 

MidS 

5.01  E+14 

+/-  1.00E+14 

0.242 

+/-  4.84E-02 

3.51  E-10 

+/- 7.01  E-11 

m 

EH 

716 

MidS 

2.00E+14 

+/-  4.00E+13 

0.112 

+/-  2.24E-02 

2.49E-10 

+/-  4.97E-11 

763 

m 

672 

EEEI 

6.12E+14 

+/-1.22E+14 

0.295 

+/-  5.90E-02 

3.85E-10 

+/-  7.71  E-11 

Table  D-3:  Emission  Indices  and  Soot  Volume  Fractions  as  a  Function  of 
Altitude,  T3  and  Fuel  -  Mid  Sulfur 
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ESP  * 

Alt* 

t3 

Fuel 

mmam 

EI(n)  unc. 

EI(m)* 

EI(m)  unc. 

PVF* 

PVF  unc. 

542 

3 

LowS 

liman 

0.00 

RSKSnSS!! 

1.548E-1 1 

+/-  3.10E-12 

545 

3 

616 

LowS 

5.113E+13 

ggEBgll 

0.04 

!fB!¥ETsfSEI 

2.236E-10 

+/-  4.47E-11 

549 

3 

672 

LowS 

6.091E+13 

gaKagaw 

0.08 

5.17E-10 

+/-  1.03E-10 

552 

3 

716 

LowS 

3.265E+13 

BBH33SH 

0.05 

3.77E-10 

+/-  7.54E-1 1 

554 

3 

794 

LowS 

1.344E+13 

+/-  2.69E+12 

0.02 

2.151  E- 1 0 

+/-  4.30E-11 

529 

9 

672 

LowS 

6.432E+13 

mbmh 

0.06 

M1MAVI 

3.515E-10 

+/-  7.03E-1 1 

532 

9 

716 

LowS 

3.924E+13 

0.05 

BBSs Jam 

3.941  E- 10 

+/-  7.88E-1 1 

534 

9 

755 

BSS9 

6.91E+12 

SMKWaMrii 1 

8.558E-1 1 

+/-  1.14E-10 

527 

9 

616 

LowS 

8.755E+13 

BBIBsBB 

0.06 

BMBaa 

1.995E-10 

+/-  3.99E-11 

498 

m 

547 

HffRJl 

4.019E+13 

0.02 

n»M»ia»a 

3.124E-1 1 

+/-  6.25E-12 

505 

12 

672 

LowS 

3.808E+13 

HMBS3SR 

0.02 

ggggga 

7.578E-1 1 

12 

716 

LowS 

5.791E+13 

BMBE3H 

0.05 

1.734E-10 

BigHin 

513 

15 

wwm 

LowS 

9.876E+13 

+/-  1.98E+13 

0.06 

1.103E-10 

+/-  2.21E-11 

514 

15 

672 

LowS 

9.384E+13 

+/-  1.88E+13 

0.06 

BMEEE1 

+/-  2.50E-11 

517 

15 

716 

LowS 

2.928E+13 

0.02 

3.605E-1 1 

+/-  7.21E-12 

* 

ESP  Escort  parameter  /  number 
Alt  Altitude 

El  (n)  Number-based  emission  index 
El  (m)  Mass-based  emission  index 
PVF  Particle  volume  fraction 
unc  uncertainty  (Is) 

Table  D-4:  Emission  Indices  and  Soot  Volume  Fractions  as  a  Function  of 
Altitude,  T3  and  Fuel  -  Low  Sulfur 
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Fuel 

men 

mm 

mm 

PVF 

|PVF  unc. 

LowS 

4.84E+13 

+/■  9.58E+12 

0.04 

+/-  7.97E-03 

1.95E-10 

+/-  4.55E-1 1 

MedS 

2.01E+14 

+/-  4.02E+13 

0.11 

+/-  2.1  IE-02 

4.08E-10 

+/-  8.16E-11 

fagga 

2.51E+14 

+/•  5.03E+13 

0.12 

+/-  2.41  E-02 

5.79E-10 

+/-  1.16E-10 

IJP8+100 

2E+14 

+/■  4.00E+13 

0.09 

+/-  1.86E-02 

3.81  E-10 

+/-  7.62E-11 

* 

ESP 

Escort  parameter  /  number 

Alt 

Altitude 

EI(n) 

Number-based  emission  index 

EI(m) 

Mass-based  emission  index 

PVF 

Particle  volume  fraction 

unc 

uncertainty  (Is) 

Table  D-i 

5:  Average  Emission  Indices  and 

Volume  Fractions  For  All  Data 


Fuel 

xbar(nm)  * 

unc 

Low  S 

63.14 

14.46 

MedS 

53.93 

10.60 

HighS 

55.86 

13.88 

JP8+100 

57.75 

7.08 

* 

xbar  (nm)  Average  mean  diameter 
unc  uncertainty  (Is) 

Table  D-6:  Average  Mean  Diameters  For  All  Data 
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m/dlogx  (pg/ cm**3) 


Aerosol  size  distr,  97829lae 
AERLOG.grf-rev,  Esc  No.  700 


x  (diam,  cm) 


Figure  D-2:  Particle  mass  distributions  showing  >  90%  of  particle  mass  falls  within 
measurement  range  from  10  nm  to  300  nm  particle  diameter 
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Aerosol  size  distr,  97826laq 
AERLOG.grf-rev,  Esc  No.  603 


Figure  D-3:  Particle  mass  distributions  showing  >  90%  of  particle  mass  falls  within 
measurement  range  from  10  nm  to  300  nm  particle  diameter 
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Figure  D-4  a-d:  El  vs  T3  for  Low  Sulfur  Case  at  altitudes  a=3km,  b=9km, 

c=12km,  d=15km 
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Figure  D-5  a-f:  El  vs  T3  for  Medium  Sulfur  Case  at  altitudes  a=3km,  b=9km,  c=l  1km, 

d=12km,  e=15km,  f=17km 
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Figure  D-6  a-d:  El  vs  T3  for  High  Sulfur  Case  at  altitudes  a=3km,  b=9km, 

c=12km,  d=15km 
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Size  Distributions  -  JP8+100 
Escort  No.  803 
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Figure  D-8:  Size  distributions  for  escort  numbers  803,  806,  811 
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Size  Distributions  -  JP8+100 
Escort  No.  814 
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Figure  D-9:  Size  distributions  for  escort  numbers  814,  817,  822 
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Size  Distributions  -  JP8+100 
Escort  No.  824 
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Figure  D-10:  Size  distributions  for  escort  numbers  824,  827,  835 
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Size  Distributions  -  JP8+100 
Escort  No.  837 
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Figure  D-l  1:  Size  distributions  for  escort  numbers  837,  840, 844 
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Figure  D-13:  Size  distributions  for  escort  numbers  868,  871, 875 
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Figure  D-14:  Size  distributions  for  escort  numbers  879,  882, 886 
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Figure  D-15:  Size  distributions  for  escort  numbers  890,  894,  899 
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Figure  D-16:  Size  distributions  for  escort  numbers  569,  572, 575 
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Figure  D-19:  Size  distributions  for  escort  numbers  599, 603, 607 
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Figure  D-20:  Size  distribution  for  escort  number  611 
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Figure  D-21:  Size  distributions  for  escort  numbers  621, 624, 627 
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Figure  D-22:  Size  distributions  for  escort  numbers  630, 633, 636 
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Figure  D-25:  Size  distributions  for  escort  numbers  659,  687, 695 
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Figure  D-26:  Size  distributions  for  escort  numbers  700, 707, 711 
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Figure  D-27:  Size  distributions  for  escort  numbers  712, 716, 718 
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Size  Distributions  -  Mid  Sulfur 
Escort  No.  720 
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Figure  D-28:  Size  distributions  for  escort  numbers  720, 763, 768 
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Size  Distributions  -  Low  Sulfur 
Escort  No.  498 
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Size  Distributions  -  Low  Sulfur 
Escort  No.  505 
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Size  Distributions  -  Low  Sulfur 
Escort  No.  508 
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Figure  D-31:  Size  distributions  for  escort  numbers  498,  505, 508 
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Size  Distributions  -  Low  Sulfur 
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Size  Distribution  -  Low  Sulfur 
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Size  Distribution  -  Low  Sulfur 
Escort  No.  549 
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Figure  D-35:  Size  distributions  for  escort  numbers  549,  552, 554 
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Figure  D-36:  Size  distributions  for  the  same  operating  conditions  recorded  at  different 

times  as  referenced  by  escort  number  sets 
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APPENDIX  E 

Chemical  Ionization  Mass  Spectrometry 

John  O.  Ballenthin,  Thomas  M.  Miller,  and  A.  A.  Viggiano 
Air  Force  Research  Laboratory 


E.l.  BACKGROUND 

The  chemical  ionization  mass  spectrometer  (CIMS)  used  in  the  jet  engine  tests  at  NASA- 
LeRC  Research  Center  was  developed  and  constructed  at  the  Air  Force  Research 
Laboratory  at  Hanscom  AFB,  Bedford,  Massachusetts.  The  instrument  has  been  flown  on 
board  NASA  T-39  and  DC-8  aircraft  to  sample  the  exhaust  gases  of  other  aircraft  in 
flights  during  the  SNIF,  SUCCESS,  and  SONEX  programs  (Ref.  E-l).  The  SNIF 
program  included  flights  behind  F-16  fighters  (using  the  F-100)  jet  engine.  The  same 
instrument,  operated  in  its  electron-bombardment  ionization  mode,  has  flown  in  the  bomb 
bay  of  a  NASA  WB-57  to  sample  the  exhaust  in  the  stratosphere  of  space  shuttle  launches 
and  Titan  missile  launches.  An  earlier  version  of  the  instrument,  also  in  its  electron- 
bombardment  mode,  has  flown  twice  on  the  space  shuttle  in  connection  with  the  “shuttle 
glow”  problem,  and  still  earlier  versions  were  flown  on  balloons,  sounding  rockets,  and 
satellites. 

Modeling  of  the  results  of  the  in  situ  aircraft  exhaust  measurements  was  complicated  by 
lack  of  knowledge  of  exactly  which  trace  gases  are  present  at  the  exhaust  plane  of  the  jet 
engine — for  example,  if  H2SO4  is  observed,  and  surely  is  produced  from  SO3  and  H2O, 
how  much  of  the  precursor  SO2  gas  is  oxidized  into  SO3  inside  the  engine  (by  O  atoms)  as 
opposed  to  oxidation  by  OH  in  the  wake?  Similarly,  how  much  NO2  is  converted  into 
HNO3  in  the  engine?  Thus,  the  NASA-LeRC  tests  are  important  to,  and  complement,  our 
other  aircraft  sampling  data. 

E.2.  CIMS  INSTRUMENT 

The  CIMS  instrument  consists  of  sampling  lines,  a  flow  reactor,  an  ion  source,  a  drying 
region,  a  sampling  orifice,  a  skimmer,  ion  lenses,  an  rf  electric  quadrupole  mass 
spectrometer,  and  a  electron  multiplier,  and  is  shown  in  Figure  E-l.  The  electronics 
package  includes  dc  and  rf  power  supplies,  pulse  counting  circuitry,  and  a  microprocessor. 
The  CIMS  instrument  weighs  about  300  lb  including  pumps  and  electronics.  The  physical 
size  is  approximately  60  x  70  x  100  cm.  A  separate  100-cfrn  mechanical  pump  was 
installed  at  NASA-LeRC  to  pump  the  flow  tube  gas  at  a  velocity  (68  m/s)  sufficient  to 
keep  the  walls  flushed,  to  reduce  backgrounds  signals,  and  to  maintain  a  low  pressure  (37 
Torr)  in  the  flow  tube.  The  37-Torr  flow  tube  pressure  was  chosen  because  the  ion  signal 
levels  maximize  in  the  20-40  Torr  range,  and  the  low  pressure  allowed  us  to  pump  engine 
effluent  through  the  sampling  capillary  at  an  acceptable  rate  over  the  entire  simulated 
altitude  range  of  the  test  cell. 
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Engine  exhaust  was  sampled  through  three  extractor  probes  on  the  lower  half  of  the  probe 
rake.  The  probes  were  located  at  5.1  cm,  10.2  cm,  and  15.2  cm  below  the  engine  axis. 
Siloxane-covered  glass-lined  stainless  steel  lines  from  the  probes  were  joined  after  1-m 
travel.  From  this  point  a  single  siloxane/glass-coated  stainless  steel  line  (0.64-cm  OD,  7.3- 
m  long),  heated  to  150  °C  to  inhibit  condensation  of  exhaust  gases,  carried  the  engine 
effluent  through  the  wall  of  the  test  cell  to  the  CIMS  instrument.  The  engine  effluent 
passed  through  a  siloxane-coated  0.32-cm  OD  stainless  steel  capillary  into  the  flow  tube  of 
the  CIMS.  The  purpose  of  the  capillary  was  to  confine  the  flow  rate  of  engine  effluent  to 
<20  slm  (std.  liters  per  min.)  to  avoid  saturating  the  CIMS  signals.  The  flow  tube  gas 
consisted  of  30  Torr  N2  flowing  continuously  at  a  rate  of  about  120  slm,  2  Torr  of  02 
from  the  ion  source  and  dryer  region,  and  typically  5  Torr  of  engine  effluent.  The  flow 
tube  (3.5-cm  OD)  was  maintained  at  150  °C  to  avoid  condensation  of  the  effluent.  A 
sampling  orifice  (100-|im  diam.)  for  the  CIMS  was  located  0.3  m  downstream  of  the  entry 
point  for  the  engine  effluent  in  the  flow  tube,  near  the  wall  of  the  flow  tube.  Directly 
opposite  the  sampling  orifice  was  a  source  of  CO3"  ions.  The  CO3  was  produced  by  a 
corona  discharge  in  a  small  cup  containing  approximately  1  atm  of  O2  gas  mixed  with 
about  0.05%  C02.  The  corona  produces  O'  ions  which  rapidly  associate  with  C02  to 
produce  CO3'  ions.  The  CO3'  ions  were  injected  into  the  flow  tube  through  a  0.13-cm 
orifice,  entrained  in  a  jet  of  source  gas  (3  slm).  The  CO3  ions  traveled  across  a  diameter 
of  the  flow  tube  to  the  sampling  orifice,  reacting  with  the  engine  effluent  along  the  way. 

The  primary  reactions  for  detection  of  S02  and  HN03  were  given  in  Section  4.3.  The 
precursor  ion  CO3'  was  chosen  for  the  present  work  because  it  does  not  react  appreciably 
with  the  major  constituents  of  the  flow  tube  gas,  N2  and  02,  but  reacts  with  SO2  and 
HNO3.  CO3'  reacts  with  S02  to  form  S03\  which  immediately  clusters  with  02  to  yield 
S05'  ion  signal  at  112  amu  (Ref.  E-2).  C03‘  forms  an  adduct  with  HNO3,  at  123  amu 
(Ref.  E-3).  One  complication  is  that  C03'  associates  with  H20  molecules,  which  are 
plentiful  in  the  combustion  exhaust.  Even  at  150  °C  the  concentration  of  CO3  (H20)  may 
be  more  intense  than  that  of  C03‘.  The  hydrated  CO3'  ions  react  with  S02  and  HNO3  to 
produce  the  same  product  ions,  but  in  some  cases  with  different  rate  constants  (Ref.  E-4). 
The  mechanism  may  be  more  complex  as  well;  for  example,  CO3'  reacting  with  SO2  forms 
SO3',  and  S03'  may  become  hydrated.  A  ligand-switching  reaction  with  02  occurs,  leading 
to  the  same  S05'  product  as  found  from  unhydrated  CO3'.  Less  important,  but  still 
significant,  are  adducts  such  as  C03'(N2)  and  C03'(C02)  and  higher-order  clusters  of  H20, 
N2,  02,  and  C02.  The  hydration  of  ions,  and  clustering  in  general,  tends  to  congest  the 
mass  spectrum  at  higher  mass  numbers.  The  detection  sensitivity  that  our  CIMS  technique 
is  capable  of  (presently  10  pptv  for  1  s  integration  time)  is  achievable  only  if  the  product 
ion  does  not  fall  at  the  same  mass  number  as  one  of  the  undesirable  cluster  species.  For 
SO2,  which  leads  to  signal  at  1 12  amu  as  described  in  Sect.  4.3,  and  HNO3,  which  leads  to 
signal  at  123  amu,  there  is  no  competing  signal. 

In  contrast,  consider  our  search  for  SO3  in  the  engine  exhaust.  S03  may  be  detected  by  the 
reaction  (Ref.  E-5): 


CO3  +  SO3  ->  SO4  +  CO2, 
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which  gives  ion  signal  at  96  amu,  the  same  mass  as  C(V(H20)2.  Since  both  SO3  and  H2O 
are  exhaust  products,  it  isn’t  possible  to  separate  contributions  at  96  amu.  A  different 
precursor  ion  such  as  (possibly)  NO3'  is  needed.  Many  species,  such  as  N2O,  simply  do 
not  react  with  CO3',  or  with  extremely  low  efficiency,  and  are  undetectable  with  this 
precursor  ion.  We  also  searched  for  H2SO4,  which  yields  ion  signal  at  97  amu  (HSO4'), 
but  found  that  it  could  all  be  accounted  for  by  HS04'  produced  in  the  flow  tube  by  SO3 
(from  SO2),  or  S05‘  interacting  with  H2O.  It  is  likely  that  neither  SO3  nor  H2SO4  can 
survive  the  8.3-m  length  of  sampling  line.  In  future  engine  tests,  we  would  like  to  place 
the  CIMS  instrument  directly  inside  the  test  cell,  with  an  orifice  or  a  very  short  sampling 
line  of  a  few  centimeters. 

Ions  entering  the  CIMS  sampling  orifice  first  passed  through  a  “drying  region”  about  0.3- 
cm  long,  where  hot,  dry  O2  gas  (3  slm)  replaced  the  flow  tube  gas  and  allowed  H2O 
molecules  to  evaporate  from  the  core  ions.  The  evaporation  is  not  complete,  but  greatly 
reduces  congestion  of  the  mass  spectrum  at  high  mass  numbers.  Between  the  sampling 
orifice  and  a  skimmer  (0.13-cm  diam.  orifice),  two  350-1/s  turbomolecular  pumps  removed 
most  of  the  gas  that  passed  through  the  sampling  orifice.  (An  oil-free  diaphragm  pump 
was  used  to  back  up  the  turbomolecular  pumps.)  An  electric  potential  of  100  V  on  the 
skimmer  tended  to  keep  ions  in  the  central  jet  of  gas  that  passed  through  the  skimmer. 
The  ion  beam  next  passed  through  electrostatic  lenses  and  entered  an  rf  quadrupole  mass 
spectrometer.  A  third  turbomolecular  pump  kept  the  quadrupole  region  at  high  vacuum 
(typically  10'5  Torr).  Ions  exiting  the  mass  spectrometer  were  detected  with  a  electron 
multiplier.  A  fourth  turbomolecular  pump  maintained  the  multiplier  region  at  about  lO"6 
Torr  vacuum.  Pulses  denoting  ion  detection  were  passed  through  shaping  and  counting 
circuitry  and  then  to  a  computer  for  data  storage  and  display. 

Capacitance  manometers  were  used  to  monitor  the  flow  tube  pressure,  the  pressure  on  the 
input  side  of  the  sampling  capillary,  and  the  differential  pressure  on  a  pitot  tube  mounted 
in  the  flow  tube.  Thermocouples  or  RTD  devices  were  used  to  monitor  temperatures  of 
the  flow  tube  and  of  the  reservoir  housing  the  HNO3  calibration  source. 

N2  purge  gas  could  be  introduced  into  the  sampling  line  inside  the  test  cell,  to  clear  engine 
effluent  out  of  the  sampling  line.  The  purge  gas  was  useful  as  a  diagnostic  to  identify  the 
source  of  certain  ion  signals — that  is,  to  show  that  a  particular  peak  in  the  mass  spectrum 
was  due  to  flow  tube  gases  instead  of  engine  effluent. 

Normally,  only  high-purity  (99.999%  pure)  O2  is  used  in  the  CIMS  ion  source,  drying 
region,  and  as  the  calibration  carrier  gas.  But  the  high-purity  O2  ordered  by  NASA-LeRC 
for  these  tests  did  not  arrive  until  the  low-sulfur  Jet  A  fuel  runs  were  completed,  meaning 
that  higher  background  signals  resulted  during  the  early  days  of  the  tests. 
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E.3.  CALIBRATION  SOURCES 

The  calibration  source  for  SO2  consisted  of  a  0.5-1  cylinder  filled  to  6  atm  of  He  gas  with 
1%  S02  and  a  flow  controller  delivering  the  mixture  to  the  flow  tube  at  the  same  port 
through  which  engine  effluent  entered  the  flow  tube.  Tests  at  the  AFRL  lab  and  at  NASA- 
LeRC  (with  the  7.3-m  long,  heated  siloxane  sampling  line)  showed  no  loss  of  SO2  gas  in 
the  sampling  line.  The  flow  controller  was  set  for  1.454  seem  (std.  cm3  per  min.)  of  the 
1%  S02  mixture  for  the  NASA-LeRC  engine  tests.  Converting  to  slm  and  putting  in  the 
1%  factor  gives  an  SO2  flow  rate  of  1 .454  x  10'5  slm. 

The  HNO3  calibration  source  consisted  of  three  HNO3  permeation  tubes  which  emitted 
approximately  5,000  ng/min  of  HNO3  each  at  60  °C.  The  HNO3  calibrant  vapor  was 
carried  from  the  60  °C  reservoir  by  400  seem  of  O2  flush  gas  passing  over  the  permeation 
tubes.  The  total  flow  rate  of  HNO3  was  approximately  15,000  ng/min  or  5.333  x  10‘6  slm. 
The  word  “approximately”  is  used  because  the  HNO3  emission  rate  depends  somewhat  on 
the  equilibration  time  of  the  permeation  tube;  a  correction  has  been  applied  to  the  data 
based  on  calibration  measurements  carried  out  at  Aerodyne  Research,  Inc.,  in  April  of 
1998.  Tests  at  NASA-LeRC  showed  about  15%  loss  of  HNO3  over  the  entire  length  of 
the  sampling  line,  as  illustrated  in  Figure  E-2.  A  caveat  is  that  we  assume  there  is  no  loss 
of  HNO3  in  the  siloxane/glass-coated,  0.32-cm  OD  stainless  steel  tubing  used  to  deliver 
the  calibration  gas  to  the  inside  of  the  test  cell.  The  He/HNCb  flow  (400  seem)  was 
maintained  in  the  delivery  line  at  all  times  to  ensure  that  the  HNO3  was  in  equilibrium  with 
surfaces  of  the  delivery  line,  so  that  there  was  no  loss  of  HNO3  up  to  the  points  of  entry 
into  the  sampling  line.  During  a  calibration  event,  a  Teflon  solenoid  valve  was  actuated  to 
divert  the  flow  into  the  sample  line.  At  other  times  the  calibrant  gas  flowed  into  the 
vacuum  chamber.  Corrections  for  HNO3  loss  in  the  sampling  lines  were  made  as  outlined 
above. 


E.4.  DATA  PROTOCOL 

The  data  protocol  included  mass  scans  from  45-127  amu  with  purge  gas  “on”  and  “off”. 
Single-mass  integrations  with  calibration  gases  “on”  and  “off’  and  purge  gas  “on”  and 
“off’  were  carried  out  for  masses  96  amu  [C03'(H20)2  and  SO4],  97  amu  (HSO4),  112 
amu  (S05‘),  and  123  amu  [CO3XHNO3)].  [Because  the  purge  gas  greatly  reduced  the 
amount  of  water  vapor  in  the  flow  tube,  thus  affecting  the  distribution  of  CO3  (H20)n 
precursor  ions  in  the  reaction  zone,  the  calibration  signal  obtained  with  purge  gas  “on”  is 
only  of  diagnostic  value.]  Next,  data  for  a  number  of  notable  masses  (46,  60,  62,  64,  78, 
80,  87,  88,  94,  96,  97,  98,  105,  107,  109,  112,  115,  116,  123,  and  125  amu)  were 
integrated  for  at  least  2.5  s,  with  the  purge  gas  “on”  and  “off.”  Ion  signal  at  most  of  these 
masses  allow  us  to  monitor  the  degree  of  H20,  N2,  and  02  clustering,  as  a  diagnostic. 
Only  the  single-mass  integrations  (over  time  periods  of  10-30  s)  were  used  to  arrive  at  the 
final  results  for  emission  indexes,  for  SO2  and  HNO3. 
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E.5.  DATA  ANALYSIS 


The  CIMS  data  analysis  for  obtaining  the  emission  index  (El)  for  a  trace  gas  will  be 
illustrated  using  S02  as  an  example.  EI(S02)  is  defined  as  the  number  of  grams  of  S02  in 
the  engine  exhaust  per  kilogram  of  fuel  burned.  We  measure  the  mixing  ratio  of  S02  in 
the  exhaust  effluent  and  use  NASA  measurements  of  combustion-related  C02  to  obtain 
EI(S02).  Specifically,  EI(S02)  =  (64/12)(0.87)(10'4)[ppbv(SO2)]  /  (%C02).  The  factors 
64  and  12  are  the  molecular  mass  of  S02  and  atomic  mass  of  C,  respectively;  the  factor 
0.87  is  the  fraction  of  C  in  the  fuel;  the  factor  104  is  the  product  of  a  factor  of  103  which 
converts  g  into  kg  and  a  factor  of  10'7  which  puts  the  ppbv(S02)  and  %C02  on  the  same 
scale;  ppbv(S02)  is  our  measured  mixing  ratio;  and  %C02  is  the  NASA-LeRC 
contribution.  EI(HN03)  is  obtained  similarly,  except  that  it  is  conventional  to  report 
EI(HN03)  in  terms  of  the  equivalent  number  of  grams  of  N02  produced  in  the  exhaust; 
hence,  the  corresponding  formula  for  EI(HN03)  contains  a  leading  factor  of  (46/12), 
where  46  is  the  molecular  mass  of  N02. 

The  CIMS  mixing  ratio  (in  ppbv)  for  S02  is  determined  from  ion  signals  [“counts(eng- 
S02)”  and  “counts(cal-S02)”]  due  to  engine-related  S02  and  calibrant  S02,  respectively, 
and  from  flow  rates  of  calibrant  S02  and  engine  effluent  [“Q(cal-SO  2)”  and  “Q(effluent)” 
as  follows: 

ppbv(eng-S02)  =  109*Q(cal-SO2)*counts(eng-SO2)  /  [counts(cal-S02)*Q(effluent)]. 

The  count  rate  for  SO5'  (due  to  S02)  was  typically  75  per  10  ms  sampling  time.  The  count 
rate  for  C03‘HN03  (due  to  HN03)  was  typically  25  per  10  ms  sampling  time.  The  results 
reported  here  were  based  on  integration  times  of  about  10  s,  so  the  total  number  of  ion 
counts  in  an  S02  measurement  was  typically  75,000,  and  for  HN03,  25,000.  S02  count 
rates  for  low-sulfur  fuel  were  smaller  than  given  above,  and  those  for  high-sulfur  fuel  were 
greater.  The  factor  of  109  in  the  formula  above  puts  the  result  in  ppbv.  The  factor  Q(cal- 
SO  2)  is  set  by  a  flow  controller  that  injects  S02  calibrant  into  the  flow  tube. 

The  flow  rate  of  engine  effluent  through  the  sampling  capillary  was  measured  both  at 
NASA-LeRC  and  in  our  laboratory,  as  a  function  of  the  pressure  across  the  capillary 
(which  ranged  from  175-1500  Torr  during  engine  tests,  depending  on  the  simulated 
altitude  and  engine  power  setting).  In  this  differential  pressure  (AP)  range,  we  found 
Q(effluent)  =  0.9434  +  0.021153  *  AP,  in  slm  for  AP  in  Torr,  at  a  temperature  of 
approximately  60  C.  During  each  test  condition,  the  pressure  on  both  sides  of  the 
capillary  were  monitored  with  capacitance  manometers  and  recorded  in  the  data  stream. 

Before  the  start  of  the  NASA-LeRC  engine  tests,  we  had  only  a  rough  idea  of  the  ion 
signal  levels  we  could  expect  from  the  engine  exhaust,  and  therefore  what  dilution  of  the 
effluent  in  the  flow  tube  would  be  required  to  avoid  saturating  the  CIMS.  Early  in  the 
engine  tests,  during  low-sulfur  fuel  runs,  our  procedures  and  ion  signals  were  still  being 
optimized  for  the  NASA-LeRC  conditions  and  these  data  carry  larger  error  bars  than  later 
data.  (The  C03'  signal  strength  was  ten  times  larger  during  the  final  week  than  during  the 
first  week.)  The  final  data  runs,  for  JP-8  fuel,  consistently  yielded  large  signal  levels  and 
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routine,  reliable  background  and  calibration  measurements,  and  hence  are  the  most 
accurate.  The  CIMS  method  is  inherently  very  accurate  as  applied  here,  if  the  flow  rates 
of  calibration  gases  and  engine  effluent  are  known,  because  the  count  rates  with  and 
without  the  calibration  gas  may  be  accurately  (<1%)  measured  by  integrating  over  10-30 
s.  The  flow  rates  of  calibration  gas  and  engine  effluent  can  be  measured  with  flow  meters 
to  within  5%.  In  practice,  however,  we  found  larger  hour-to-hour  variations  in  the  results, 
as  great  as  20%,  indicating  a  systematic  error  associated  with  the  sampling  and/or 
calibration  system  (but  not  with  the  CIMS  itself).  The  systematic  error  (for  example,  a 
leak  in  the  sampling  line  or  in  the  calibration  line)  cannot  be  traced  now  that  the  tests  have 
ended.  Considerable  effort  has  gone  into  tracking  down  the  source  of  the  problem, 
unsuccessfully.  However,  it  is  clear  for  several  reasons  that  future  tests  should  be 
conducted  with  the  CIMS  mounted  inside  the  test  cell,  with  short  sampling  lines.  Because 
of  the  systematic  error,  the  present  results  are  only  accurate  to  ±30%.  Averages  for 
EI(S02)  over  many  days  of  data  are  estimated  accurate  within  15%  for  JP-8  fuel,  20%  for 
medium-sulfur  Jet  A  fuel,  and  25%  for  low-sulfur  Jet  A  fuel.  There  is  greater  variation  in 
EI(HNC>3)  because  of  sampling  line  losses  and  associated  time  constants. 


E.6.  RESULTS 

E.6.1.  Nitric  Acid  Production  In  the  Engine 

NO2  and  OH,  produced  in  the  combustion  process,  combine  to  produce  HNO3.  In  this 
section  we  will  examine  the  efficiency  of  that  process.  While  we  do  not  have  a 
measurement  of  OH  concentration,  the  NASA-LeRC  measurements  of  NO  and  NOx 
concentration  allow  us  to  deduce  the  N02  concentration.  The  emission  indexes  for  NO, 
NOx,  and  HNO3  are,  by  convention,  all  expressed  in  terms  of  grams  of  NO2  per  kilograms 
of  fuel.  Therefore,  we  may  calculate  EI(NC>2)  =  EI(NOx)  -  EI(NO),  and  compare  this 
value  to  our  EI(HN03)  measured  under  the  same  engine  conditions.  Figure  E-3  shows 
EI(N02)  plotted  versus  combustor  inlet  temperature  (T3)  for  all  fuels  and  all  altitudes  for 
which  we  have  corresponding  EI(HN03)  data  (i.e.,  not  all  of  the  NASA-LeRC  NO2  data 
are  represented).  Much  of  the  dispersion  in  the  data  shown  in  Figure  E-3  is  due  to  an 
altitude  dependence  in  EI(N02);  there  is  little  dependence  on  fuel  type. 

In  Figure  E-4  are  shown  the  fraction  of  NO2  conversion  into  HN03,  calculated  from 
EI(HN03)  /  [EI(N02)  +  EI(HN03)].  The  dispersion  in  the  data  now  includes  the  scatter  in 
our  measurements  of  EI(HN03).  The  most  striking  feature  of  Figure  E-4  is  the  rapid 
decline  in  the  production  of  HNO3  as  the  combustor  inlet  temperature  increases.  The  N02 
— »  HNO3  conversion  rate  decreases  a  factor  of  5  between  600-700  °K,  and  drops  at  least  a 
decade  between  700-800  °K.  This  decline  is  not  due  to  a  corresponding  change  in  the 
concentration  of  N02;  indeed,  EI(N02)  increases  by  perhaps  a  factor  of  3  in  the  600-800 
°K  range  (Figure  E-3).  We  attribute  the  decrease  in  the  N02  -4  HN03  conversion  rate  to 
the  thermal  instability  of  HNO3. 

In  Figures  E-5  and  E-6  we  examine  the  altitude  dependence  in  these  results,  for  one  fuel 
type  (#2).  In  Figure  E-5,  the  N02  -4  HNO3  conversion  rate  is  seen  to  depend  on  altitude. 
We  shall  show  in  the  following  sections  that  there  is  no  significant  altitude  dependence  in 
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the  EI(HN03)  measurements.  Figure  E-6  shows  that  the  altitude  dependence  in  the 
calculation  of  the  N02  — >  HN03  conversion  rate  is  a  consequence  of  division  by 
EI(N02)],  which  is  altitude  dependent.  If  EI(HN03)  does  not  track  EI(N02)  precisely,  it 
indicates  that  some  other  aspects  of  the  conversion  reaction  are  changing  with  altitude, 
e.g.,  EI(OH),  the  temperature  profile  in  the  combustor,  and  the  residence  time  in  the 
combustor. 

Figures  E-7  through  E-13  give  the  N02  -*  HN03  conversion  rate  for  each  altitude  for  all 
fuel  types  for  which  we  obtained  data  on  EI(HN03).  A  linear  fit  to  each  data  set  is 
indicated  on  the  figures  to  facilitate  intercomparisons.  N02  -*  HN03  conversion  rates  for 
each  fuel  type  will  be  given  in  the  following  sections.  Overall  averages  are  8%  at  500  °K, 
6%  at  600  °K,  1.5%  at  700  °K,  and  <0.1%  at  800  °K. 

E.6.2.  JP-8  FUEL 

We  will  first  focus  on  the  JP-8  results  because  our  procedures  and  signal  strengths  were 
well  optimized  during  the  final  week  of  the  test  series.  The  emission  indexes  for  S02 
production,  EI(S02),  are  presented  in  Figs  E-14  through  E-20,  as  a  function  of  combustor 
inlet  temperature  (T3).  EI(S02)  is  judged  to  be  independent  of  combustor  inlet 
temperature  at  all  altitudes,  within  the  scatter  in  the  data.  Complementary  plots  vs 
combustor  inlet  pressure  (P3)  are  not  given  because  there  is  a  similar  lack  of  dependence 
on  P3.  The  value  of  EI(S02)  averaged  over  combustor  inlet  temperature  is  somewhat 
different  at  different  altitudes;  it  is  our  opinion  that  these  variations  are  artifacts  of  the 
measurement  system,  as  outlined  in  the  previous  section.  Thus,  it  is  best  to  consider  the 
average  (over  altitude  and  combustor  inlet  temperature)  of  all  measurements  of  EI(S02) 
for  JP-8  fuel,  <EI(S02)>  =  0.573  ±  0.086  g/kg. 

Fuel  sulfur  analyses  were  obtained  from  Wright -Patterson  AFB  (336  ppmm  for  JP-8).  If 
all  the  sulfur  (32  g/mol)  were  converted  to  S02  (64  g/mol),  the  fuel  would  yield  an  S02 
emission  index  0.672  g/kg.  Assuming  the  fuel  sulfur  analyses  are  correct,  we  find  an 
overall  average  of  (85  ±  15)%  for  the  amount  of  sulfur  appearing  as  gas  phase  SO2.  The 
remaining  sulfur,  if  any,  is  most  likely  in  the  form  of  S03  or  is  incorporated  into 
particulates. 

The  HN03  emission  indexes  (expressed,  by  convention,  in  terms  of  grams  of  N02 
produced)  are  given  in  Figures  E-21  through  E-27.  The  amount  of  HN03  produced  is 
found  to  be  independent  of  the  sulfur  level  in  the  fuel,  and  decreases  with  combustor  inlet 
temperature  as  outlined  in  the  previous  section.  For  JP-8,  we  find  EI(HN03)  =  0.20 
g(N02)/kg  fuel  at  500  °K.  This  figure  drops  to  0.17  g(N02)/kg  at  600  °K,  and  decreases 
still  further  to  about  0.05  g(N02)/kg  at  700  °K. 

* 

E.6.3.  MEDIUM-SULFUR  JET  A  FUEL 

As  with  the  other  fuels,  there  is  no  significant  dependence  of  EI(S02)  on  combustor  inlet 
temperature  for  the  medium-sulfur  (152  ppmm)  Jet  A  fuel.  The  results  are  given  in 
Figures  E-14  through  E-20.  EI(S02)  is  similarly  independent  of  combustor  inlet  pressure 
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(P3)  within  the  scatter  in  the  data.  The  overall  average  EI(S02)  for  the  medium-sulfur  Jet 
A  fuel  is  0.337  ±  0.067  g/kg.  The  fractional  uncertainty  for  the  medium-sulfur  fuel  results 
is  somewhat  larger  than  for  the  JP-8  because  signal  levels  were  lower  and  the  S02 
calibration  is  less  certain  (the  calibration  mixture  is  no  longer  available  to  verify).  The  fuel 
sulfur  analysis  from  WPAFB  implies  that  the  maximum  EI(S02)  is  0.304  g/kg.  Assuming 
this  latter  figure  is  accurate,  the  correct  EI(S02)  must  be  in  the  range  0.270-0.304  g/kg. 
Thus,  our  best  estimate  is  that  88-100  percent  of  the  fuel  sulfur  appears  in  the  form  of  S02 
in  the  engine. 

EI(HN03)  behaves  as  with  the  other  fuels:  EI(HN03)  =  0.12  ±  0.05  g(N02)  per  kg  fuel  at 
500  °K.  This  number  drops  by  10%  at  600  °K,  then  drops  a  further  factor  of  10  by  750 
°K.  The  data  are  presented  in  Figures  E-21  through  E-27. 

E.6.4.  HIGH-SULFUR  JET  A  FUEL 

Again,  EI(S02)  does  not  appear  to  vary  with  engine  combustor  inlet  temperature  (T3) 
within  the  scatter  in  the  data,  here  for  the  high-sulfur  (1113  ppmm)  Jet  A  fuel.  The  results 
are  given  in  Figures  E-14  through  E-20.  EI(S02)  is  likewise  independent  of  combustor 
inlet  pressure  (P3)  within  the  scatter  in  the  data.  The  overall  average  EI(S02)  for  the 
medium-sulfur  Jet  A  fuel  is  2.49  ±  0.62  g/kg.  The  fractional  uncertainty  is  still  larger  for 
this  fuel  because,  at  that  point  in  the  test  series,  we  were  just  settling  on  routine 
operation — with  flow  tube  oxygen  quality,  flow  tube  heating,  and  calibrant  gas  delivery 
system  all  being  modified  between  data  runs.  The  fuel  sulfur  analysis  from  WPAFB  implies 
that  the  maximum  EI(S02)  is  2.226  g/kg.  If  so,  the  correct  EI(S02)  must  be  in  the  range 
1.87-2.23  g/kg.  Thus,  our  best  estimate  is  that  84-100  percent  of  the  fuel  sulfur  appears  in 
the  form  of  S02  in  the  engine  exhaust. 


The  EI(HN03)  data  for  the  high-sulfur  fuel  show  greater  statistical  variation  than  for  the 
fuels  discussed  above,  but  the  gross  behavior  is  the  same.  EI(HN03)  =  0.15  ±  0.05 
g(N02)  per  kg  fuel  at  500  K.  This  number  drops  to  0.13  at  600  °K,  then  drops  a  decade 
by  750  K.  The  data  are  presented  in  Figures  E-21  through  E-27. 

E.6.5.  LOW-SULFUR  JET  A  FUELS 

The  pre-test  analyses  of  the  two  batches  of  low-sulfur  Jet  A  fuel  at  WPAFB  yielded  “<1 
ppmm”  sulfur  content.  Our  data  for  the  mixing  ratios  of  S02  in  the  engine  exhaust  for  all 
fuels  indicate  that  the  first  batch  of  low-sulfur  fuel  was  actually  3.5  ppmm  sulfur,  and  the 
second  batch  was  7.4  ppmm  sulfur.  These  results  represent  extrapolations  of  our 
measurements  of  mixing  ratios  versus  the  WPAFB  analyses  for  the  two  medium-sulfur 
fuels  (one  Jet  A  and  one  JP-8)  and  for  the  high-sulfur  fuel.  The  sulfur  concentrations  given 
here,  3.5  and  7.4  ppmm,  respectively,  are  technically  lower  limits  because  in  deriving  these 
values  we  are  assuming  that  the  percentage  of  fuel  sulfur  appearing  as  S02  in  the  exhaust 
was  the  same  for  the  low-sulfur  fuels  as  measured  (85-100%)  for  the  higher-sulfur  fuels. 
(It  would  seem  unreasonable,  to  us,  for  the  low-sulfur  fuels  to  have  a  lower  conversion 
rate.) 
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Results  for  EI(SC>2)  are  given  in  Figures  E-14  through  E-20.  There  is  no  dependence  on 
combustor  inlet  temperature  (or  pressure,  P3)  within  the  scatter  in  the  data.  Because  these 
data  are  from  the  first  week  of  the  test  series,  when  signal  strengths  were  low,  other 
parameters  had  not  been  optimized,  and  the  data  protocol  was  still  being  developed,  the 
uncertainties  in  even  the  average  values  for  EI(S02)  are  relatively  large  (25%).  The  first 
batch  of  low-sulfur  fuel  yielded  (EI(SC>2))  =  0.007  ±  0.002  g/kg.  The  second  batch  of  low- 
sulfur  fuel  yielded  <EI(S02)>  =  0.015  ±  0.004  g/kg. 

EfrHNCb)  follows  the  pattern  seen  with  the  other  fuels,  though  the  data  for  the  low-sulfur 
fuels  have  the  greatest  scatter,  again  due  to  the  lower  signal  levels  in  the  earlier  data  runs. 
We  find  EI(HN03)  =  0.12  ±  0.06  g(NC>2)  per  kg  fuel  at  500  °K,  dropping  to  about  0.10  at 
600  °K,  and  then  declining  more  than  a  decade  by  750  °K  (see  Figures  E-21  through  E- 
27). 


E.7.  CONCLUSIONS 

We  monitored  SO2  and  HNO3  emissions  in  the  F-100  engine  exhaust  with  the  CIMS 
instrument,  and  measured  mixing  ratios  for  these  compounds.  Comparing  to  NASA-LeRC 
measurements  of  C02  production  by  the  engine,  the  mixing  ratios  were  given  as  emission 
indexes  in  Figures  E-14  through  E-20,  for  the  various  fuels  (Jet  A  with  low-sulfur, 
medium-sulfur,  and  high-sulfur,  and  JP-8)  and  altitudes  ranging  from  3-17  km,  as 
functions  of  engine  combustor  inlet  temperature. 

The  emission  indexes  for  SO2  are  essentially  independent  of  altitude  and  of  combustor 
inlet  temperature  and  pressure.  Combining  our  measurements  of  SO2  mixing  ratios, 
NASA-LeRC  data  on  CO2  emissions,  and  sulfur  analyses  of  the  fuels  carried  out  at 
Wright-Patterson  AFB,  we  find  that  85-100%  of  the  fuel  sulfur  appears  as  gaseous  S02  in 
the  engine  exhaust.  Within  our  experimental  uncertainty,  this  conversion  level  is 
independent  of  fuel  type,  altitude,  and  power  setting.  The  remaining  fuel  sulfur,  if  any, 
most  likely  appears  as  gas  phase  SO3  or  is  incorporated  into  particulates. 

The  amount  of  HNO3  produced  is  independent  of  the  sulfur  level  in  the  fuel.  The  HNO3 
level  decreases  considerably  with  engine  temperature,  presumably  because  HNO3  is  not 
stable  at  high  temperatures.  Within  the  uncertainties  in  the  data,  there  is  no  difference  in 
HNO3  production  for  the  various  fuels  studied.  Using  NASA-LeRC  measurements  of 
EI(NOx)  and  EI(NO),  and  the  AFRL  measurements  of  EI(HN03),  we  can  calculate  the 
NO2  — >  HNO3  conversion  rate  in  the  engine.  The  NO2  — »  HNO3  conversion  rate  is  about 
8%  at  500  K,  6%  at  600  K,  1.5%  at  700  K,  and  drops  at  least  a  decade  by  800  °K.  (These 
figures  are  averages  over  all  altitudes  and  all  fuels  studied.)  We  attribute  this  decrease  in 
the  NO2  — >  HNO3  conversion  rate  to  the  thermal  instability  of  HNO3.  There  is  a 
systematic  altitude  dependence  in  this  conversion  rate.  A  detailed  examination  of  the  data 
(Figures  E-3  through  E-13)  leads  us  to  conclude  that  the  altitude  dependence  is  due  either 
to  variations  in  EI(OH) — not  measured  in  this  test  series — or  due  to  altitude-dependent 
changes  in  combustor  temperature  profile  and/or  residence  time  affecting  HNO?  formation 
from  OH  and  NO2.  The  steepness  of  the  EI(HN03)  curve  versus  T3  will  be  a  critical  test 
for  engine  combustion  models. 
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Figure  E-l.  AFRL  chemical  ionization  mass  spectrometer. 
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Figure  E-2.  Ion  signal  due  to  HNO3  injected  sequentially  at  three  points  along  the 
sampling  line.  The  data  indicate  that  there  is  a  15%  loss  of  HNO3  over  the 
full  8.3-m  length  of  the  sampling  line. 
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Figure  E-3.  EI(N02)  fromEI(NOx)  -EI(NO),  for  all  fuels  and  altitudes. 
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Figure  E-4.  Fraction  of  N02  converted  into  HNO3,  for  all  fuels  and  altitudes. 
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Figure  E-5.  Fraction  of  NO2  converted  into  HNO3  for  fuel  2  at  different  altitudes. 
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Figure  E-6.  EI(N02)  for  fuel  2  deduced  from  NASA-LeRC  data 
for  EI(NOx)  and  EI(NO). 
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Figure  E-7.  Fraction  of  NO2  converted  into  HNO3  at  3  km  altitude. 


Figure  E-8.  Fraction  of  NO2  converted  into  HNO3  at  9. 1  km  altitude. 
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E-9.  Fraction  of  NO2  converted  into  HNO3  at  10.7  km  altitude. 


Figure  E-10.  Fraction  of  NO2  converted  into  HNO3  at  12.2  km  altitude. 
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Figure  E-l  1.  Fraction  of  NO2  converted  into  HNO3  at  13.7  km  altitude. 


Figure  E-12.  Fraction  of  NO2  converted  into  HNO3  at  15.2  km  altitude. 
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Figure  E-13.  Fraction  of  NO2  converted  into  HNO3  at  16.8  km  altitude. 
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Figure  E-14. 
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Figure  E-21.  EI(HNC>3)  at  3  km  altitude  (simulated  sea  level). 
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Figure  E-24.  EI(HN03)  at  12.2  km  altitude. 
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Figure  E-26.  EI(HNO;s)  at  15.2  km  altitude. 
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APPENDIX  F 

Infrared  Tunable  Diode  Laser  System 

J.  Wormhoudt,  T.  Berkoff  and  R.  C.  Miake-Lye 
Center  for  Aero-Thermodynamics 
Aerodyne  Research,  Inc.,  Billerica,  MA  01821-3976 


F.l.  TDL  INSTRUMENT 

Tunable  lead  salt  diode  lasers  have  a  long  history  of  spectroscopic  applications',  and  have  recently 
been  applied  with  considerable  success  to  sensitive  detection  of  atmospheric  trace  gases.2  Our 
system  had  been  applied  to  combustion  gas  flows  in  four  previous  field  tests,  first  in  a  combustor 
simulator  at  NASA  Lewis  Research  Center  (LeRC)3,  then  in  two  engine  tests  at  the  Air  Force 
Arnold  Engineering  Development  Center  (AEDC)45,  and  finally  in  a  1996  engine  test  in  the 
Propulsion  Systems  Laboratory  test  cell  PSL-4  which  also  supported  the  test  reported  here. 

The  laser  diodes  are  housed  in  a  liquid  nitrogen  dewar,  with  their  nominal  operating  temperature 
controlled  by  a  small  resistance  heater.  Coarse  frequency  tuning  is  achieved  by  changing  the  diode 
temperature  and  thereby  its  refractive  index  and  cavity  mode  frequencies.  Fine  tuning  is  done  by 
scanning  the  injection  current,  which  also  has  the  effect  of  slightly  changing  the  diode  temperature. 
This  current  scanning  can  be  done  rapidly,  allowing  the  averaging  of  many  spectra  whose  spectral 
features  are  unperturbed  by  low-frequency  fluctuations  in  the  transmission.  In  the  work  reported 
here,  the  number  of  spectral  points  acquired  was  500  and  the  scanning  rate  was  300  kHz  (spectral 
points  per  second)  so  that  600  scans  of  each  diode  were  accumulated  each  second. 

A  schematic  drawing  of  the  laser  optical  table  is  given  in  Figure  F-l.  The  cryogenic  dewar  houses 
four  laser  diodes  behind  each  of  two  infrared  transmitting  windows,  with  translation  of  the  dewar  or 
the  collection  optics  serving  to  select  which  diode  emits  into  the  collecting  optics  and  so  through  the 
entire  optical  train  until  laser  light  is  focused  onto  a  detector  located  in  the  same  dewar.  The  optical 
system  collects  light  from  two  laser  diodes,  sending  two  coincident,  collimated  beams  to  a  either  a 
retroreflector  or  multiple  pass  cell  in  the  altitude  test  chamber,  then  collecting  the  return  beams  and 
focusing  them  onto  the  same  detector  element.  Measurement  time  scales  allow  the  use  of  a  single 
detector,  with  discrimination  between  the  two  lasers  achieved  by  alternating  laser-on  (frequency 
scanning)  and  laser-off  periods.  This  design,  allowing  the  simultaneous  operation  of  two  laser  diodes, 
originated  in  a  desire  to  collect  data  as  efficiently  as  possible  during  engine  tests.  It  can  also  be  a  key 
requirement  for  the  line-ratio  emission  index  measurement  technique  whose  investigation  was  the 
focus  of  the  previous  test  and  several  of  the  first  days  of  the  current  test.  In  those  cases,  the  two 
diodes  were  used  to  measure  NO  and  CO2  lines.  The  main  focus  of  the  present  test  series  was  the 
measurement  of  the  sulfur  oxide  species  SO2  and  SO3.  One  diode  was  used  to  measure  each  species, 
with  water  lines  present  in  each  spectrum  serving  as  reference  lines  for  emission  index  calculations. 

We  begin  a  discussion  of  the  details  of  the  optical  layout  by  tracing  the  paths  of  the  main  or 
diagnostic  beams.  Infrared  light  from  the  laser  diodes  is  collected  by  reflective  microscope 
objectives  (15x)  and  focused  onto  200  pm  pinholes,  which  define  input  apertures.  The  pinholes  are 
used  only  during  alignment,  so  they  are  mounted  on  removable  kinematically  indexed  bases.  The 
microscope  objectives  are  mounted  on  X-Y-Z  translators  to  allow  accurate  focusing  into  the  fixed 
apertures.  Past  the  input  aperture,  curved  mirrors  create  approximately  collimated  beams.  The  two 
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laser  beams  are  then  combined  at  the  main  beamsplitter,  and  the  co-aligned  beams  pass  to  a  pair  of 
flat  mirrors  and  onto  a  large  gimbal-mounted  steering  mirror  which  directs  the  beam  off  the  table 
and  into  the  engine  test  cell.  An  important  development  for  this  test  was  the  replacement  of  the 
former  main  beamsplitter,  a  relatively  thick  (3  mm)  sheet  of  ZnSe,  with  a  very  thin  pellicle  which 
eliminated  one  source  of  interference  fringes  which  can  be  a  major  noise  source  in  our  spectra.  The 
return  beam  is  collected  by  the  large  steering  mirror  and  a  curved  mirror  of  similar  diameter,  passed 
through  several  more  flat  mirrors  and  finally  focused  onto  the  main  beam  detector. 

The  remaining  combined  beam  associated  with  the  main  beamsplitter  goes  into  the  reference  optical 
path,  where  it  may  be  sent  through  a  gas  calibration  cell  and/or  passed  through  a  grating 
monochromator,  then  focused  onto  the  reference  detector.  Some  calibration  cells  contain  a  sufficient 
pressure  of  the  gases  of  interest  that  their  strong  absorption  features  fully  absorb  the  laser  light.  This 
makes  it  easy  to  recognize  that  the  laser  is  operating  in  the  right  wavelength  region,  and  to  check  for 
the  presence  of  other  laser  modes.  These  absorption  lines  are  also  used  to  set  the  frequency  scale  of 
the  laser  scan.  In  this  case,  we  used  cells  containing  NO  and  N02  as  well  as  N20,  to  locate  lines  in 
the  C02  region,  and  a  cell  of  methane  to  determine  frequency  scales  in  the  two  sulfur  oxide  spectral 
regions.  The  monochromator  is  used  to  provide  a  coarse  measurement  of  the  laser  frequency  during 
setup  of  the  instrument,  and  to  determine  whether  the  diode  is  emitting  at  different  wavelengths 
(different  modes).  Since  no  monochromator  is  available  in  the  main  beam  path,  diode  operating 
conditions  must  be  chosen  which  allow  single  mode  operation. 

There  is  a  parallel  visible  optical  system  for  alignment  and  setup.  A  red  HeNe  "trace"  laser  beam 
passes  through  the  dichroic  beamsplitter  downstream  of  each  microscope  objective,  and  is  coaligned 
with  each  infrared  beam.  Coalignment  is  guaranteed  by  focusing  the  beam  through  the  input 
aperture.  The  trace  beam  is  an  indispensable  aid  for  alignment  of  the  optical  system.  In  addition,  the 
trace  beam  is  used  for  accurate  calibration  of  the  monochromator,  via  higher  order  diffraction.  The 
fourth  port  of  the  dichroic  beamsplitter  can  be  used  to  observe  the  laser  diode  during  alignment.  An 
eyepiece  mounted  at  the  position  conjugate  to  the  pinhole  forms  an  effective  150X  microscope. 

We  found  we  had  to  make  one  addition  to  the  optical  system  due  to  an  effect  which  was  presumably 
present  in  the  previous  PSL  test  but  not  noticed  then  due  to  the  more  forgiving  nature  of  the 
retroreflector  optical  system.  The  effect  was  the  motion  of  the  test  cell  and  the  optics  inside  with 
respect  to  the  laser  table  outside  the  cell  which  was  firmly  attached  to  the  concrete  floor.  This 
motion  first  occurred  as  the  test  cell  was  brought  to  a  simulated  altitude,  and  seemed  roughly 
proportional  to  the  test  cell  static  pressure.  The  maximum  displacement  was  about  0.6  cm.  Initially, 
we  compensated  for  this  by  moving  the  entire  laser  table  whenever  there  were  major  changes  in  the 
test  cell  pressure.  Later,  we  installed  a  translation  mount  on  the  steering  mirror  on  the  laser  table. 
This  made  recovery  of  a  good  transmitted  signal  relatively  straightforward. 

The  temperature  and  current  of  each  laser  are  controlled  by  a  Laser  Photonics  controller.  All  inputs 
to  this  controller  were  made  via  a  program  running  on  the  data  acquisition  and  analysis  computer.  In 
a  typical  experiment,  the  laser  temperature  is  held  constant  while  the  current  through  the  laser  is 
modulated  with  a  computer  generated  sawtooth  to  sweep  the  output  frequency  across  the  infrared 
transition.  The  sawtooth  was  generated  in  250  discrete  steps  for  each  laser  by  a  digital  to  analog 
converter  board.  The  total  of  500  display  points  in  the  two  laser  spectra  is  swept  at  a  rate  of  300 
kHz  (display  points  per  second).  During  approximately  ten  percent  of  the  duty  cycle  the  laser 
current  is  dropped  below  the  lasing  threshold  to  provide  a  precise  measurement  of  zero  light 
intensity. 
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The  data  acquisition  techniques,  involving  a  Scientific  Solutions  Lab  Master  analog  to  digital 
converter  board  and  direct  memory  access  transfer  to  the  extended  memory  of  the  data  acquisition 
computer,  remain  as  described  previously.3  Although  the  TDL  instrument  could  in  principle  produce 
real-time  records  of  gas  concentrations,  this  would  require  input  of  probe  rake  data  (or  assumptions 
about  properties  along  the  line  of  sight  from  other  sources).  Our  data  acquisition  therefore  consisted 
of  the  recording  of  a  large  number  of  spectral  files  covering  segments  of  both  the  S02  and  S03 
regions  (or  the  NO  and  C02  regions).  These  were  recorded  as  averages  over  periods  of  10  to  600  s. 


F.2.  DIAGNOSTIC  IMPLEMENTATION  IN  ENGINE  TEST  CELL 

The  optical  system  used  to  propagate  the  laser  beam  through  the  engine  exhaust  in  the  altitude  test 
cell  was  substantially  modified  from  the  one  used  in  the  previous  test  in  PSL-4.  The  TDL  instrument 
was  again  placed  on  a  table  directly  in  front  of  an  optical  access  port  into  the  test  cell.  However, 
this  table  rested  on  the  concrete  floor  of  the  test  cell  building,  where  the  earlier  table  had  been  on  the 
mezzanine  above  it.  Therefore,  the  TDL  beams  were  directed  into  the  test  cell  in  the  horizontal 
plane,  rather  than  downward  (at  approximately  a  45  degree  angle)  as  in  the  1996  test. 

A  change  was  necessitated  by  the  larger  size  of  the  engine  in  the  current  test,  and  the  ground  level 
mounting  seemed  to  offer  fewer  problems  with  vibration  as  well  as  a  simpler  optical  layout.  A  0.95 
cm  thick  CaF2  window  with  a  1°  wedge  (to  suppress  interference  fringes)  was  mounted  on  the  access 
port  using  a  retaining  plate  with  a  10  cm  diameter  clear  area.  The  access  port  had  a  direct  view  of 
the  engine  exit  plane  location,  so  that  the  laser  beam  could  pass  through  the  exhaust,  strike  a  comer 
cube  retroreflector  mounted  on  the  test  cell  floor,  and  return,  horizontally  displaced  by  about  2  cm, 
to  the  test  cell  access  port  window  and  the  steering  mirror.  The  HeNe  trace  beam  could  be  followed 
through  the  entire  optical  path,  which  allowed  precise  aiming  of  each  steering  mirror. 

However,  most  of  our  observations  in  this  test  series  were  made  not  with  the  retroreflector  but  with 
a  multipass  cell  which  caused  the  laser  beams  to  pass  through  the  exhaust  14  times  before  being 
returned  to  the  optical  table.  The  multipass  mirrors  were  15  cm  diameter,  91.4  cm  focal  length 
mirrors,  with  the  mirror  nearest  the  laser  having  a  1  cm  hole  in  its  center.  They  were  contained 
inside  heavy  aluminum  boxes  equipped  with  nitrogen  purge  lines.  The  gaseous  nitrogen  purge  was 
needed  not  only  as  a  protection  against  recirculating  exhaust  but  also  to  remove  as  much 
atmospheric  water  vapor  and  methane  from  the  path  as  possible.  The  path  between  the  access  port 
and  the  back  of  the  near  box  was  contained  inside  a  10  cm  diameter  aluminum  tube,  also  purged 
with  gaseous  nitrogen.  Stainless  steel  tubes,  30  cm  in  length  and  15  cm  in  diameter,  were  attached 
to  the  front  of  each  protective  box.  This  left  a  122  cm  gap  between  the  two  tubes.  In  the  center  of 
this  gap  was  the  exhaust  flow,  exiting  from  a  nozzle  with  a  diameter  of  about  65  cm.  The  separation 
between  the  outer  two  of  four  probe  rake  thermocouples  which  register  the  full  core  exhaust 
temperature  was  30.48  cm,  while  the  pair  of  thermocouples  separated  by  45.72  cm  read 
temperatures  more  than  half  way  to  ambient  values.  The  separation  between  the  fronts  of  the 
protective  boxes,  where  thermocouples  were  mounted,  was  182  cm. 
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F.3.  EXAMPLES  OF  OBSERVED  SPECTRA 


Figure  F-2  presents  an  example  of  raw  diode  intensity  spectra  for  the  TDL  diagnostic  while  the  laser 
beam  was  passing  through  the  exhaust.  The  vertical  axis  is  simply  intensity  at  the  detector 
expressed  in  mV.  The  horizontal  axis  is  in  display  points,  or  data  acquisition  points  over  the  laser 
scan.  The  left  hand  side  of  the  scan  is  the  intensity  of  the  S02  laser.  After  this  laser  is  turned  off  at 
250  display  points  and  the  intensity  briefly  drops  to  near  the  zero  level,  the  S03  laser  is  turned  on, 
and  its  intensity  is  recorded  in  the  second  half  of  the  scan.  The  strongest  absorption  feature  seen  in 
the  S02  laser  scan,  at  display  point  210,  is  a  water  absorption  line  at  1332.757  cm'1. 

In  the  S03  side  of  the  laser  scan,  a  diode  mode  break  (abrupt  transition  in  the  laser  frequency)  occurs 
at  315  display  points.  The  region  between  250  and  315  display  points,  containing  some  very  strong 
water  lines,  was  not  used  in  our  analysis.  Only  the  region  between  315  and  490  display  points  was 
analyzed,  and  even  there  the  absorption  features  seen  in  the  figure  are  water  lines.  The  smaller  and 
more  numerous  S03  lines  lie  under  and  between  the  water  lines.  An  analysis  was  carried  out  on  the 
highest  altitude  spectra  for  the  high  sulfur  loading,  where  S03  peaks  should  be  most  clearly 
distinguishable.  By  averaging  absorbances  in  S03  peak  and  valley  regions,  we  were  able  to 
determine  that  our  detection  limit  for  S03  was  less  than  2  ppmv,  and  that  the  exhaust  concentrations 
of  S03,  even  under  these  most  favorable  conditions,  were  below  that  level. 


F.4.  DATA  ANALYSIS  PROCEDURES  AND  ERROR  LIMIT  ESTIMATES 
F.4.1.  Data  Analysis  Using  the  TDL  Data  Acquisition  Program 

Our  technique  in  the  data  analysis  reported  here  is  to  measure  relative  column  densities  of  water  and 
S02  and  then  use  this  ratio  and  the  known  emission  index  for  water,  derived  from  the  fuel  formula, 
to  yield  a  direct  measurement  of  S02  emission  index  without  passing  through  the  steps  of 
quantifying  absolute  molecular  mixing  fractions  and  mass  fluxes.  To  carry  out  the  relative  column 
density  determination,  the  absorption  features  of  water  and  S02  are  quantified  by  least  squares  fits, 
using  the  TDL  data  acquisition  program  in  data  analysis  mode.  This  is  the  routine  which  allows  real¬ 
time  fitting  and  concentration  measurement,  here  applied  to  stored  spectra.  Lines  are  fit  to  a  set  of 
Voigt  line  profiles  using  a  nonlinear  least  squares  routine  based  on  the  Levenberg-Marquardt 
approach.6 

A  fit  to  the  entire  region  seen  on  the  left  half  of  Figure  F-2  is  dominated  by  the  strong  water  lines  at  the 
right  hand  side  of  the  S02  scan.  To  quantify  the  water  relative  column  density,  we  actually  used  a  more 
restricted  region,  from  108  to  248  display  points.  The  S02  determination  was  made  using  a  fit  to  one  of 
several  even  more  restricted  regions,  detailed  below.  The  baseline,  varying  from  unit  transmission  due  to 
diode  laser  power  fluctuations  and  residual  etalon  fringes,  was  represented  by  a  fourth  order  polynomial 
over  each  region  in  turn.  The  line  parameters  used  in  the  least  squares  fits  are  taken  from  the  HITRAN 
compilation7  or  from  more  recent  measurements  taken  from  the  literature  or  communicated  to  us  by  the 
investigators  prior  to  publication.  The  absorption  line  properties  are  evaluated  at  a  single  representative 
temperature  and  pressure  derived  from  averages  of  probe  rake  values.  The  TDL  analysis  also  uses  an 
input  path  length  (here  taken  as  500  cm)  to  derive  an  absolute  density  along  this  uniform  absorption  path. 
This  path  length,  14  times  the  distance  through  the  hot  region  of  the  plume  as  estimated  from  probe  rake 
data,  will  disappear  when  the  column  density  ratios  are  formed,  but  it  does  in  fact  lead  to  absolute  mixing 
fractions  which  agree  well  with  sampling  values. 
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Figure  F-3  shows  the  entire  spectral  region  which  was  analyzed.  In  this  figure  the  solid  line  is  the 
observed  trace  as  converted  to  a  transmission  scale  by  the  TDL  data  analysis  program,  through  its  fit 
of  a  polynomial  baseline.  The  frequency  scale  is  set  using  reference  spectra  of  the  methane 
reference  cell,  taken  in  between  exhaust  spectra.  The  dashed  line  is  the  actual  least  squares  fit  by  the 
TDL  data  analysis  program,  evaluated  at  the  same  data  acquisition  points.  This  spectrum,  for  high 
sulfur  loading,  a  50  kft  simulated  altitude  and  a  T3  value  of  629  °K,  has  lines  which  are  fairly  well 
separated  from  each  other.  As  the  altitude  drops  and  the  pressure  rises,  the  lines  become  broader 
and  overlap  each  other  more,  and  it  becomes  more  difficult  to  make  an  accurate  estimate  of  the 
position  of  the  unit  transmission  baseline.  As  described  above  a  fit  of  the  region  from  1332.45  cm'1 
to  1332.85  cm'1,  including  the  strong  lines  on  the  right  hand  side,  is  used  to  determine  the  water 
concentration. 

Figure  F-4  shows  the  spectral  region  containing  the  spectral  features  used  to  determine  the  S02 
concentration.  The  circles  show  the  positions  of  actual  data  acquisition  points.  The  two  strongest 
absorption  lines  in  this  spectrum,  below  1332.1  cm’1  and  above  1332.4  cm’1,  are  due  to  methane  in 
the  test  cell  air  in  the  regions  between  the  exhaust  and  the  purged  boxes  containing  the  optics.  The 
smaller  peak  near  1332.3  cm'1  is  due  to  water,  which  also  has  a  few  much  smaller  lines  scattered 
through  this  region.  Most  of  the  remaining  peaks,  in  particular  those  near  1332.2  cm'1  and  between 
1332.35  and  1332.4  cm'1,  are  due  to  S02.  This  is  seen  in  Figure  F-5,  which  shows  model  spectra  for 
the  individual  species  S02  and  H20  as  well  as  their  sum  and  the  observed  spectrum  of  Figure  F-4. 
When  the  S02  features  are  strong  enough,  it  is  possible  to  fit  spectra  of  all  three  species  to  the  entire 
spectral  region  in  Figure  F-5  to  derive  three  species  concentrations.  When  S02  features  are  weaker, 
it  is  necessary  to  concentrate  on  smaller  spectral  ranges  in  order  to  maximize  the  accuracy  of  the 
derived  S02  concentration.  Therefore,  for  most  observations  of  the  highest  sulfur  loading,  the 
spectral  region  which  was  fit  to  determine  S02  concentrations  was  from  1331.9  to  1332.45  cm1,  (29 
to  108  display  points),  while  for  the  lower  sulfur  spectra,  the  most  common  range  used  was  from 
about  1332.15  to  1332.35  cm'1  (51  to  88  display  points).  This  latter  range  is  seen  to  cover  a  spectral 
region  including  the  single  water  line  just  below  1332.3  cm1,  and  the  S02  feature  (composed  of 
several  lines)  around  1332.2  cm'1. 


F.4.2.  Derivation  of  Emission  Index  Values 

As  noted  above,  the  absolute  concentration  values  resulting  from  TDL  data  analysis  are  not  entirely 
measured  quantities,  but  instead  incorporate  several  approximations.  These  include  the  assumptions 
of  an  arbitrary  absorption  path  length  and  of  uniform  pressures  and  temperatures  along  the  line  of 
sight.  These  pressures  and  temperatures  are  derived  from  averages  of  probe  rake  values  for  those 
sampling  points  fully  in  the  exhaust  core  flow.  For  analysis  of  infrared  data,  the  array  of  Mach 
numbers  from  probe  measurements  and  an  assumed  specific  heat  ratio  for  the  exhaust  flow  are  used 
to  convert  probe  rake  observed  total  temperatures  to  static  temperatures,  while  a  similar  conversion 
is  applied  to  measured  total  pressures  to  yield  an  array  of  static  pressures.  Finally,  ratios  of  total 
temperatures  in  the  exhaust  and  the  test  cell  free  stream  are  used  to  estimate  mixing  ratio  profiles, 
which  are  then  used  as  weights  in  averages  of  the  static  pressure  and  temperature  arrays  derived 
from  probe  rake  data  to  give  single  pressures  and  temperatures  for  input  to  the  TDL  program.  After 
carrying  out  a  few  calculations  using  a  multi-cylinder  model  whose  properties  are  determined  by 
individual  probe  values,  we  estimate  the  error  levels  associated  with  the  assumption  of  uniform 
average  properties  as  being  on  the  order  of  a  few  per  cent,  when  mixing  ratio  weighted  averages  are 
used. 
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Emission  index  is  defined  as  the  weight  in  grams  of  an  exhaust  component  per  kilogram  of  fuel 
burned.  Thus,  for  the  major  exhaust  species  C02  and  water,  which  are  produced  solely  and 
quantitatively  from  the  hydrogen  and  carbon  in  the  fuel  (combining  with  oxygen  from  the  air), 
emission  indices  are  determined  once  the  chemical  composition  of  the  fuel  is  known.  For  the 
primary  fuel,  involved  in  sulfur  loading  variations,  four  analyses  were  made,  with  an  average  value 
for  water  emission  index  of  1240  (and  a  standard  deviation  of  less  than  5).  For  the  alternate  fuel 
(JP-8  +100)  we  used  the  first  value  we  received  of  1269,  even  though  a  second  analysis  yielded  a 
value  of  1253.  Therefore,  for  the  primary  fuel  case,  the  formula  used  is 

1240  (64/18)  ([S0J/[H20]) 

where  64/18  converts  from  mass  of  H20  to  mass  of  S02.  Since  we  know  the  weight  fraction  of 
sulfur  in  each  fuel,  we  can  calculate  maximum  theoretical  S02  emission  indices  for  the  case  when  all 
fuel  sulfur  is  present  in  the  exhaust  as  S02.  Since  the  molecular  weight  of  S02  is  twice  that  of  atomic 
sulfur,  these  maximum  theoretical  S02  emission  indices  are  simply  twice  the  weight  fraction  of 
sulfur  in  the  fuel  (again,  expressed  as  grams  per  kilogram).  For  the  three  fuels  for  which  we  could 
measure  S02,  these  maximum  emission  indices  were  2.226  for  high  sulfur,  0.672  for  JP-8,  and  0.304 
for  the  mixture  of  high  sulfur  and  low  sulfur  fuels. 


F.4.3.  Measured  S02  Emission  Index  Values 

We  will  present  S02  emission  index  values  in  a  series  of  six  tables.  More  than  one  table  is  required 
because  we  want  to  examine  not  only  the  absolute  values  of  S02  emission  index  for  each  observation 
condition,  but  also  their  ratios  to  theoretical  maximum  values  (fractional  conversion  of  fuel  sulfur  to 
S02)  and  their  averages  over  one  or  two  of  the  three  observation  parameters  (sulfur  loading,  power 
setting,  and  altitude).  The  error  limits  in  the  tables  have  very  different  meanings,  as  well,  again 
because  there  are  several  issues  we  want  to  examine. 

One  difference  in  the  error  limits  in  the  various  tables  stems  from  the  conventional  division  of  error 
limit  components  into  two  types,  random  and  systematic  errors.  Estimates  of  random  errors  are 
typically  derived  from  the  statistics  of  multiple  measurements,  while  estimates  of  systematic  errors 
involve  bounding  the  effects  of  as  many  potential  sources  of  bias  or  uncertainty  as  possible. 
Especially  when  small  numbers  of  measurements  are  taken,  estimates  of  random  errors  may  vary 
widely  from  one  measurement  condition  to  another,  while  estimates  of  fractional  systematic  errors 
could  either  be  made  specific  to  each  measurement  condition  or  (our  choice  in  the  estimates  described 
in  the  following  section)  could  be  assumed  to  apply  equally  to  all  measurements.  The  error  limits  in 
Tables  F-2,  F-4  and  F-6  are  total  uncertainty  limits  which  incorporate  estimates  of  both  types  of  error. 
Our  work  to  date  in  making  estimates  of  systematic  errors  is  detailed  in  the  following  section.  The 
result  documented  there  was  a  decision  to  adopt  a  constant  fractional  error  limit  of  0.16  as  an 
estimate  of  systematic  errors.  This  value  is  added  in  quadrature  to  relative  standard  deviations  in 
multiple  measurements  (derived  from  the  absolute  standard  deviations  in  Table  F-l  by  dividing  by  the 
maximum  theoretical  emission  index  values)  to  form  the  error  limits  quoted  in  Table  F-2. 

The  other  difference  in  the  error  limits  is  related  to  the  question  of  what  one  wishes  to  define  as  a 
single  measurement  of  a  quantity,  such  as  S02  conversion  fraction.  The  error  limits  in  Table  F-l  are 
simply  the  standard  deviations  in  arrays  of  emission  index  measurements  determined  from  individual 
spectra.  In  the  cases  where  only  one  spectrum  was  obtained,  no  error  limit  is  listed.  Such  error  limits 
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provide  useful  information  about  the  reproducibility  of  some  measurements,  and  the  variability  of 
others,  at  the  spectrum-to-spectrum  level.  Another  point  of  view,  however,  is  that  the  measurement 
of  S02  conversion  fraction  is  not  complete  until  several  spectra  have  been  obtained  and  an  average 
value  has  been  formed.  In  that  view,  the  averages  in  Table  F-l  become  the  primary  measured  data, 
and  it  is  then  of  interest  to  see  what  the  random  variations  are  in  those  values.  In  Tables  F-3  and  F- 
5,  then,  the  error  limits  are  standard  deviations  representing  the  distributions  of  emission  index 
values  in  Table  F-l. 

We  can  now  return  to  Table  F-l,  which  presents  S02  emission  index  values  for  all  sulfur  loadings, 
power  settings,  and  altitudes  for  which  we  could  make  determinations  with  reasonable  accuracy  (we 
will  comment  below  on  the  remaining  cases,  those  involving  the  lowest  sulfur  loading  and  the  lowest 
simulated  altitude).  For  convenience,  we  also  divide  by  the  theoretical  maximum  emission  index 
values  and  present,  in  Table  F-2,  the  fraction  of  fuel  sulfur  in  the  exhaust  in  the  form  of  S02.  It  can  be 
seen  that  unity  (all  fuel  sulfur  existing  as  S02  in  the  exhaust)  lies  within  or  close  to  the  range  of 
variation  for  all  cases. 

The  data  in  Tables  F-l  and  F-2  serve  as  the  inputs  to  averages  used  to  generate  the  remaining  tables. 
In  these  tables,  averages  over  one  or  more  of  the  three  variables  are  carried  out  in  order  to  make 
clearer  any  possible  trends  with  changes  in  the  remaining  variables.  For  instance,  Tables  F-3  and  F-4 
can  be  used  to  consider  the  possibility  of  a  systematic  trend  with  simulated  altitude.  In  fact,  they 
clearly  demonstrate  that  within  our  uncertainty  limits,  we  do  not  observe  any  trend  with  altitude. 

In  Tables  F-5  and  F-6,  all  data  points  with  the  same  simulated  altitude  have  been  averaged,  allowing 
consideration  of  trends  with  power  setting  and  sulfur  loading.  It  can  be  seen  that  in  most  cases,  the 
emission  indices  at  the  lowest  and  highest  power  settings  are  smaller  than  those  in  the  middle.  It  is 
entirely  possible  that  this  is  an  experimental  artifact,  due  to  two  trends  with  engine  power  setting 
which  make  the  measurements  at  the  lowest  and  highest  values  the  most  difficult  to  perform.  At  idle, 
the  species  concentrations  and  line  absorbances  are  smallest  and  least  easily  measured  against  baseline 
fluctuations.  At  maximum  power,  greater  fluctuations  in  laser  beam  position  on  the  detection  optics 
reduced  the  intensity  of  transmitted  laser  light,  again  lowering  the  accuracy  of  the  measurement. 

The  final  trend  to  be  considered,  that  with  sulfur  loading,  is  most  easily  considered  using  the  final 
column  in  Table  F-6.  Although  the  S02  conversion  fraction  seems  to  increase  substantially  with 
decreasing  fuel  sulfur  concentration,  it  can  also  be  seen  that  it  can  be  taken  to  be  constant  within  the 
estimated  error  limits.  For  the  three  fuel  types,  the  absolute  values  of  S02  emission  indices  averaged 
over  all  engine  conditions  and  altitudes  are  1.64  !  0.18,  0.55  !  0.11  and  0.27  !  0.05.  The  error  limits 
here  reflect  only  the  standard  deviations  in  individual  data  points.  Dividing  by  theoretical  upper  limits 
to  the  emission  indices,  they  yield  fractional  values  of  0.08,  0.16  and  0.16  respectively.  The  overall 
average  of  all  the  individual  data  points  of  Table  2  yields  the  sulfur  conversion  fraction  in  the  bottom 
right  comer  of  Table  F-6  of  0.82  !  0.22.  This  value  and  its  error  limits  (now  reflecting  systematic  as 
well  as  random  errors)  encompass  most  of  the  individual  values  in  Table  F-2. 

The  data  sets  in  which  we  have  the  most  confidence  are  those  at  the  highest  altitudes.  It  can  be  seen 
that  the  sulfur  conversion  fractions  are  lowest  for  the  higher  altitude  data  sets.  However,  the  fact  that 
average  values  of  S02  fraction  are  higher  at  the  lowest  altitudes  is  most  likely  not  significant  and 
simply  an  artifact  of  the  greater  difficulty  of  the  measurements.  Given  the  uncertainty  estimates,  in  no 
case  can  we  say  that  some  fuel  sulfur  is  definitely  not  in  the  form  of  S02  in  the  exhaust.  Instead,  the 
most  we  can  say  about  the  possibility  of  other  sulfur  species  in  the  exhaust  is  that  they  could  exist,  but 
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only  at  levels  of  30  per  cent  of  total  fuel  sulfur  at  the  very  most.  In  the  specific  case  of  S03,  this  upper 
limit  is  consistent  with  our  upper  limit  of  less  than  2  ppmv:  at  high  altitudes,  this  concentration  level 
of  S02  converts  to  an  emission  index  increment  of  about  0.25.  This  increment  is  1 1  per  cent  of  the 
theoretical  maximum  for  the  high  sulfur  loading  fuel.  In  other  words,  both  our  S02  and  S03 
measurements  are  consistent  with  the  statements  that  at  least  the  substantial  majority,  and  quite 
possibly  the  overwhelming  majority,  of  fuel  sulfur  exists  in  the  form  of  S02. 

Finally  we  turn  to  those  components  of  the  data  set  which  are  not  reported  in  the  tables.  The  lowest 
simulated  altitude,  leading  to  exhaust  static  pressure  values  in  the  500  to  1000  Torr  range,  results  in 
absorption  lines  so  overlapped  that  an  accurate  determination  of  the  baseline  position  has  so  far  not 
been  possible.  The  existence  of  S02  absorption  features  is  discernible,  and  plausible  baseline  values 
can  be  seen  to  be  consistent  with  S02  levels  at  the  theoretical  maximum.  However,  though  we  can 
show  that  these  low  altitude  data  subsets  are  consistent  with  values  derived  at  higher  altitudes,  we 
cannot  use  them  to  derive  independent  measurements. 

The  low  sulfur  fuel  was  expected  to  be  well  below  our  S02  detection  limit.  However,  we  did  perform 
an  analysis  of  a  subset  of  this  data  which  led  to  an  estimate  of  the  minimum  detectable  S02  loading  of 
less  than  1  ppmv.  This  minimum  detectable  concentration  should  be  smaller  than  the  random 
fluctuations  in  larger  concentrations  to  be  discussed  in  the  following  section,  and  indeed  it  will  be 
seen  that  the  random  error  component  of  the  total  error  estimate  is  in  the  range  of  1  to  3  ppmv  for 
high  sulfur  loading  data  points,  and  around  1  ppmv  for  the  mixed  sulfur  data  points. 


F.4.4.  Estimates  of  the  Errors  in  Measured  S02  Emission  Index  Values 

We  almost  always  took  at  least  two  spectra  for  each  engine  power  setting.  The  standard  deviations 
in  emission  indices,  when  available,  are  taken  to  represent  the  random  error  component  of  our 
uncertainty  levels.  We  observed,  as  we  have  in  past  tests,  that  reproducibility  between  parameters 
obtained  from  successive  spectra  is  often  very  good.  However,  the  spectral  data  analysis  for  SO2  is 
a  more  demanding  problem  that  was  the  case  of  NO,  for  example,  and  we  do  see  more  spectrum-to- 
spectrum  variation  in  derived  concentrations  than  in  past  tests.  In  several  cases  observations  were 
made  of  the  same  nominal  engine  condition  on  two  different  days.  We  observed  that  while  the  engine 
parameters  were  very  accurately  reset,  as  seen  in  very  similar  probe  rake  data,  the  variation  between 
our  observed  concentration  values  for  different  days  was  larger  than  the  variations  in  a  single  day. 
Because  the  range  of  standard  deviations  was  large  and  because  systematic  errors  turned  out  to  be  the 
dominant  contributors  to  overall  error  bounds  in  most  cases,  we  did  not  include  a  random  error 
component  in  the  estimates  when  only  one  spectrum  was  used.  The  average  standard  deviations  in 
Table  F-l  are,  for  high  sulfur  data  points,  0.16  or  a  fractional  standard  deviation  of  0.07,  0.05 
absolute  or  0.08  relative  for  JP-8,  and  0.07  absolute  or  0.23  relative  for  mixed  sulfur  fuel.  These 
random  error  limits  could  also  be  assigned  to  the  two  values  in  Table  F-l  which  represent  single 
spectra.  In  Table  F-5,  in  which  standard  deviations  are  now  formed  from  the  distributions  in  average 
values  of  Table  F-l,  a  similar  procedure  could  be  followed  for  the  highest  power  setting  values,  thus 
also  raising  the  total  error  limits  in  this  column  of  Table  F-6. 

The  estimation  of  possible  systematic  errors  was  carried  out  by  identifying  those  input  parameters  to 
the  data  analysis  program  which  have  significant  uncertainties,  perturbing  them  by  that  uncertainty 
limit,  and  carrying  out  a  determination  of  S02  emission  index.  For  some  parameters,  the  issues 
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involved  with  choosing  an  uncertainty  limit  and  determining  the  sensitivity  of  the  emission  index  to 
such  a  variation  are  completely  straightforward,  while  for  others  they  are  more  complex.  An 
example  of  the  latter  situation  is  the  temperature  along  the  absorption  path.  The  data  analysis 
program  assumes  a  single  temperature  along  the  entire  path.  In  fact,  the  temperature  varies 
substantially  through  the  exhaust.  Because  of  this,  a  number  of  systematic  errors  can  result  from  the 
assumption  of  a  single  temperature.  As  mentioned  above,  we  carried  out  a  few  calculations  in  which 

*  we  determined  the  S02  emission  index  using  from  a  multi-cylinder  model  which  incorporated  all  the 
measurements  from  the  probe  rake.  We  then  found  that  this  emission  index  was  close  to  that 
derived  from  a  single-cylinder,  or  uniform  path,  model  whose  temperature  was  derived  from  the 

'  array  of  probe  rake  temperatures  by  weighted  averaging  using  a  mixing  ratio  derived  from  ratios  of 

total  temperatures.  The  disagreement  between  detailed  and  single-temperature  models  is  not  the  only 
potential  systematic  error  we  wish  to  estimate.  For  example,  the  finite  spatial  resolution  of  the 
probe  rake  could  lead  to  uncertainties  in  temperature  over  significant  segments  of  the  total  path. 

As  an  easily  available  measure  of  the  uncertainties  associated  with  the  temperature  input  to  the 
analysis  program,  we  examined  the  differences  in  temperatures  derived  from  two  different  averaging 
procedures.  One  was  the  mixing  ratio  weighted  technique  mentioned  above,  while  the  other  was  a 
simple  linear  average  of  the  four  innermost  temperatures.  On  the  average,  these  two  values  differed 
by  about  40  °K,  a  value  we  adopted  as  the  uncertainty  in  the  temperature  input  to  the  analysis 
program.  The  sensitivity  of  S02  and  H20  concentrations  and  the  resulting  emission  index  turns  out 
to  be  quite  variable,  since  this  parameter  can  influence  the  spectral  fit  in  a  variety  of  ways. 
Averaging  sensitivity  calculations  carried  out  using  a  half  dozen  spectra  covering  a  range  of 
altitudes,  power  settings  and  sulfur  concentrations  had  the  result  that  a  40  °K  perturbation  in  the  path 
temperature  had  a  9  per  cent  effect  on  the  S02  emission  index. 

Two  points  should  be  made  here.  The  first  is  that  uncertainties  in  individual  probe  rake 
measurements  undoubtedly  are  much  smaller  than  the  above  !40  °K,  which  is  intended  to  represent 
an  entire  complex  of  uncertainties  associated  with  the  analysis  of  absorption  over  a  nonuniform  path. 
In  an  earlier  discussion  of  rake  data  uncertainty5,  a  typical  fractional  uncertainty  in  a  measured  probe 
[total]  temperature  was  said  to  be  0.0075,  resulting  in  a  an  absolute  uncertainty  in  this  temperature 
regime  of  the  order  of  !5  °K.  The  discussion  goes  on  to  say  that  the  maximum  fractional  uncertainty 
in  any  of  the  derived  thermodynamic  parameters  due  to  propagation  of  measurement  uncertainties 
through  the  data  reduction  algorithms  is  about  0.02,  and  even  when  calibration  uncertainties  are 
included  the  total  fractional  uncertainty  in  any  rake  aerodynamic  property  is  conservatively  estimated 
to  be  no  more  than  0.05. 

The  second  point  to  be  made  is  that  even  a  9  per  cent  sensitivity  to  a  40  °K  change  in  temperature 

y  could  be  substantially  reduced  if  S02  and  H20  spectral  features  could  be  used  whose  strength 

variation  with  temperature  was  better  matched  than  those  used  here.  In  this  first  measurement  of 
S02  in  an  aircraft  engine  exhaust  by  infrared  absorption,  we  were  constrained  by  the  laser  diodes 

*  available  to  us  at  the  time  to  the  spectral  regions  reported  here.  In  preparation  for  future 
observations,  we  plan  to  investigate  other  spectral  regions.  There  is  the  possibility  that  the  estimated 
sensitivity  to  uncertainties  in  the  characterization  of  the  temperature  along  the  absorption  could  be 
reduced  to  a  negligible  contribution  to  the  overall  error  limit  estimate,  through  the  use  of  different 
spectral  features. 

If  the  same  type  of  estimate  is  made  for  an  uncertainty  in  the  pressure  input  to  the  analysis  program, 
using  the  differences  between  the  two  types  of  averaging,  we  find  an  average  fractional  uncertainty 
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of  about  0.05.  If  we  examine  the  variation  of  static  pressure  values  across  the  flow,  we  find  average 
fluctuations  closer  to  a  fraction  of  0.10,  even  though  the  true  static  pressure  should  be  close  to 
constant.  However,  even  this  larger  uncertainty  level  leads  to  an  uncertainty  in  S02  emission  index 
of  only  5  per  cent,  which  when  added  in  quadrature  with  several  larger  contributions  has  only  a 
small  effect  on  the  total  uncertainty  estimate.  An  even  more  negligible  source  of  error  turned  out  to 
be  uncertainty  in  the  diode  tuning  rate.  Calibration  spectra  were  taken  every  few  minutes,  and  the 
tuning  rate  derived  from  the  closest  one  is  normally  used  to  analyze  each  spectrum.  Taking  the 
tuning  rate  from  the  next-closest  reference  spectrum  resulted  in  changes  in  S02  emission  index 
which  were  at  most  1  per  cent,  and  usually  much  less. 

Our  fitting  procedure  in  analyzing  this  data  set  involved  keeping  the  model  for  the  zero  absorption 
baseline  fixed  at  a  fourth  order  polynomial,  then  varying  the  beginning  and  ending  points  of  the 
fitted  region  and  the  exact  position  (in  display  points)  of  a  reference  point  for  the  absolute  frequency 
scale  (in  the  S02  region,  the  water  peak  at  1332.29  cm'1  is  used)  until  a  good  visual  fit  is  obtained  to 
the  S02  peak  at  1332.2  cm'1.  The  height  of  this  absorption  peak  is  typically  only  about  twice  the  size 
of  the  minimum  absorption  on  either  side  of  it.  In  other  words,  in  this  spectral  region  we  cannot 
directly  determine  the  position  of  the  baseline.  In  analyses  of  earlier  data  sets  in  which  the  baseline 
position  was  better  determined  (and  so  the  above  procedure  of  varying  input  parameters  until  a  good 
visual  fit  was  obtained  was  not  used),  we  estimated  the  error  due  to  this  uncertainty  by  varying  the 
polynomial  order.  In  the  data  set  of  Reference  5,  this  was  done  for  all  NO  spectra,  with  typical 
variations  being  in  the  range  of  a  few  per  cent.  For  the  current  S02  data  set,  we  only  performed  this 
set  of  analyses  on  one  representative  spectrum,  finding  about  a  5  per  cent  effect  for  reasonable 
variations  of  polynomial  order.  Again,  this  level  of  sensitivity  makes  it  a  relatively  small 
contributor  to  the  overall  uncertainty  estimate. 

Finally,  we  consider  uncertainties  in  spectroscopic  parameters,  whose  effects  on  the  emission  index 
turn  out  to  be  straightforward,  giving  very  similar  sensitivities  in  all  spectra  studied.  Our  S02  line 
parameters  were  taken  not  from  the  HITRAN  listing7  but  from  a  more  recent  reanalysis8.  This  work 
found  that  HITRAN  line  strength  values  were  often  10  to  25  per  cent  stronger  than  their  values,  due 
to  spectroscopic  perturbations  ignored  in  HITRAN  and  to  differences  in  the  band  strengths  used. 
We  could  take  a  substantial  fraction  of  these  differences  as  characterizing  the  remaining  uncertainty 
in  S02  line  strengths.  However,  differences  in  our  spectral  region  turn  out  to  be  only  about  3  per 
cent,  so  we  adopt  this  value,  leading  directly  to  a  3  per  cent  change  in  S02  emission  index  and  thus 
only  a  modest  contribution  to  the  total  systematic  error  limit.  We  assume  that  uncertainties  in  water 
line  strengths  are  even  smaller. 

The  water  line  width  values  we  used  are  based  on  recent  measurements  communicated  to  us 
privately  by  R.  A.  Toth  of  the  Jet  Propulsion  Laboratory,  taking  into  account  calculations  of 
temperature  dependence  performed  for  us  by  R.  R.  Gamache  of  the  Center  for  Atmospheric 
Research  of  the  University  of  Massachusetts,  Lowell.  Perhaps  the  most  useful  estimate  of 
uncertainty  in  these  parameters  is  the  observation  that  even  when  experimental  errors  are  assessed 
by  the  groups  that  made  them  as  being  only  a  few  per  cent,  comparisons  between  different 
measurements  suggest  that  true  error  limits  (including  systematic  errors  which  are  harder  to  assess) 
are  closer  to  10  per  cent.9  Using  this  value  leads  consistently  to  variations  of  about  6  per  cent  in  S02 
emission  index.  We  should  point  out  that  the  10  per  cent  real  error  limit  referred  to  room 
temperature,  measured  values,  while  we  must  use  values  extrapolated  to  the  500  K  range,  if  anything 
adding  to  their  uncertainty. 
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The  S02  room  temperature  line  width  values  we  used  were  all  0.10  cm'1,  reduced  from  the  HITRAN 
value  of  0.12  to  achieve  better  fits  with  our  observed  spectra.  Independently,  Sumpf,  Schoene  and 
Kronfeldt  measured  a  number  of  nearby  lines  (though  none  of  the  lines  we  used)  and  reported  an 
average  line  width  for  our  band  of  0.10  ±  0.01 1  cm'1.  Again,  because  we  must  take  into  account  the 
further  uncertainty  of  the  line  width  value  at  elevated  temperatures,  we  use  this  entire  error  limit  in 
our  sensitivity  calculations,  arriving  at  an  average  9  per  cent  change  in  S02  emission  index  for  SO 
line  widths  perturbed  by  1 1  per  cent  from  the  nominal  value.  2 

In  Table  F-7,  we  summarize  the  above  discussion  of  possible  systematic  errors,  and  show  that  the 
overall  estimate  for  fractional  systematic  errors  is  0.19.  This  is  comparable,  and  if  anything 
somewhat  smaller,  than  estimates  for  previous  data  sets.  Furthermore,  it  is  clear  that  with  additional 
work  it  can  be  reduced.  Ratioing  S02  and  water  lines  with  more  comparable  temperature 
dependences  has  already  been  mentioned  as  a  way  of  lowering  the  sensitivity  to  probe  rake 
temperatures.  Additional  measurements  of  spectroscopic  parameters  could  also  help.  In  addition, 
the  random  error  component  measured  by  standard  deviations  in  multiple  measurements  could  be 
reduced  through  better  understanding  of  the  fluctuations  in  the  baseline,  both  those  due  to  turbulence 
which  are  very  strong  for  short  averaging  times,  and  those  due  to  interference  fringes  which  remain 
even  after  averaging  for  many  minutes.  This  leads  us  to  the  conclusion  that  measurements  with  total 
error  limits  under  10  per  cent  are  well  within  the  realm  of  possibility.  Even  the  total  error  limits 
shown  in  the  tables,  in  the  range  of  17  to  25  per  cent,  allow  useful  conclusions  to  be  drawn  from  the 
current  data  set. 


F.5.  SUMMARY  AND  CONCLUSIONS 

In  this  section  we  wish  to  address,  quantitatively,  the  key  issues  involved  in  understanding  our  data 
set  of  S02  emission  indices  or  S02  conversion  fractions.  One  set  of  issues  involves  the  fact  that  S02 
conversion  fraction  essentially  has  both  an  upper  and  a  lower  bound.  The  upper  bound  is  definite: 
the  sulfur  in  the  exhaust  cannot  exceed  the  sulfur  in  the  fuel.  The  lower  bound  comes  from  our 
upper  limit  on  S03  concentration  and  the  assumption  that  other  sulfur  species  are  similarly  small 
components  of  the  total  sulfur  in  the  exhaust.  Our  goal  in  this  area  is  to  demonstrate  that  the  absolute 
values  in  our  data  set  and  our  estimates  of  error  limits  are  consistent,  so  that  none  of  these  bounds 
are  exceeded  by  more  than  our  estimated  uncertainties.  A  second  set  of  issues  involves  the 
assumption  that  systematic  errors  can  be  reduced  as  the  measurement  technique  is  further  developed, 
so  that  it  is  of  interest  to  examine  the  precision,  as  opposed  to  the  accuracy,  of  the  S02  conversion 
fractions  determined  by  our  data  set.  We  will  discuss  each  area  in  turn. 

Our  consistency  with  the  upper  limit  constraint  on  S02  conversion  fraction  has  already  been 
discussed.  Table  F-2  shows  a  few  values  which  exceed  1,  but  none  by  more  than  our  error  estimates. 

*  The  consistency  with  a  lower  limit  constraint  is  also  easily  demonstrated.  We  recall  that  our 

examination  of  the  SO3  detection  limit  in  the  high  sulfur,  highest  altitude  spectra  yielded  an  upper 
limit  of  2  ppmv  which  was  equivalent  to  an  S02  emission  index  of  about  0.25,  or  an  S02  conversion 
fraction  of  0.11.  It  happens  that  these  high  altitude,  high  sulfur  spectra  also  give  the  lowest  values  of 
S02  conversion  fraction.  This  can  be  seen  in  Table  F-4,  which  also  lists  the  appropriate  total  error 
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limits.  If  we  take  the  55  kft  data  point,  the  sum  of  S02  conversion  fraction  plus  total  uncertainty  limit 
plus  upper  limit  for  SO3  is 

0.72  +  0.19  +  0.11  =  1.02 

so  that  our  S02  and  S03  measurements  and  error  limits  are  indeed  consistent. 

We  can  estimate  the  ultimate  precision  inherent  in  a  data  set  of  the  size  discussed  here  by  computing 
the  standard  deviation  of  the  mean  of  S02  conversion  fraction.  This  is  not  the  standard  deviation  of 
the  distribution,  which  should  remain  roughly  constant  as  more  samples  are  considered,  but  is 
obtained  from  it  by  dividing  by  the  square  root  of  the  number  of  samples,  so  that  as  the  number  o 
samples  increases,  the  precision  of  the  mean  value  increases  (even  though,  due  to  systematic  errors, 
the  measured  value  may  not  be  the  true  value).  If  we  consider  each  spectrum  to  be  an  independent, 
equivalent  measurement  of  S02  conversion  fraction,  then  we  have  130  samples  and  while  the 
standard  deviation  in  S02  conversion  fraction  for  this  entire  array  is  0.23,  the  standard  deviation  in 
the  mean  is  only  0.020.  If  we  choose  to  consider  the  three  fuels  as  three  different  data  sets,  thereare 
still  30  to  more  than  50  samples  in  each,  and  standard  deviations  in  the  mean  are  0.14  for  hig 
sulfur,  0.026  for  JP-8  and  0.036  for  the  mixed  sulfur  fuel. 

These  levels  of  precision  are  of  considerable  interest  because,  as  discussed  above,  we  believe  it  is 
possible  with  additional  work  to  reduce  the  sytematic  errors  present  in  our  first  S02  data  set.  ine 
importance  of  better  understanding  of  these  systematic  errors  can  be  understood  by  considering  the 
average  values  for  sulfur  conversion  fraction  for  the  three  fuels,  shown  in  the  last  column  in  Tab  e 
F-6.  These  values,  0.74,  0.82  and  0.89  going  from  highest  to  lowest  sulfur,  clearly  differ  by  much 
more  than  the  precisions  in  mean  values  quoted  above.  There  are,  of  course,  two  extreme 
possibilities:  there  is  a  real  change  in  sulfur  conversion  fraction  with  fuel  sulfur  loading,  or  there  are 
systematic  effects  on  our  measurements  resulting  in  three  different  values  when  in  fact  there  is  only 
one.  There  could  also  be  both  a  real  dependence  on  sulfur  loading  and  systematic  errors  as  well. 
Our  exercise  in  estimating  possible  systematic  errors  turns  out  not  to  allow  a  choice  between  these 
possibilities,  since  it  only  shows  that  systematic  errors  could  be  as  large  as  the  differences  between 
sulfur  loading  data  sets.  Another  possible  reason  for  rejecting  a  systematic  error  would  be  that  the 
trend  does  not  logically  follow  from  a  change  in  the  spectra  with  sulfur  loading.  However,  it  is  not 
impossible  that  as  the  S02  peaks  become  comparable  to  incompletely  characterized  background 
peaks  of  water  and  methane,  more  of  the  intensity  in  some  background  peaks  is  taken  by  the  least 
squares  fit  to  be  due  to  S02,  thus  giving  rise,  on  the  average,  to  the  apparent  trend  quoted  above. 

In  summary,  tunable  diode  laser  measurements  in  the  S02  and  S03  spectral  regions  were  carried  out 
for  the  first  time  in  an  aircraft  engine  exhaust  and  give  the  same  result:  fuel  sulfur  is  conve  e 
largely  to  S02  at  the  exhaust  exit  plane,  with  S02  forming  at  least  70  per  cent  of  the  total  sulfur 
species  under  all  conditions,  with  typical  fractions  in  the  80  to  90  per  cent  range. 
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Table  F-l.  S02  Emission  Index  Values  for  Individual  Engine  Power  Settings,  Fuel  Sulfur  Loadings 
and  Simulated  Altitudes,  with  Standard  Deviations  in  Multiple  Measurements 


Power 

Setting, 

Altitude,  m 

T2JK 

16,764 

15,240 

13.716 

12.192 

10.668 

9144 

HighS 

500-538 

1.4110.14 

1.3110.06 

1.6010.11 

572 

1.85 

616-630 

1.8110.12 

1.8710.25 

1.5810.05 

672-686 

1.76±0.22 

1.6210.20 

1.5810.04 

2.0210.27 

1.6310.13 

714-720 

1.4410.09 

1.5410.17 

1.5110.21 

1.7610.32 

1.7110.13 

783 

1.4910.12 

JP-8 

505-506 

0.4810.09 

0.3410.16 

583-619 

0.5410.04 

0.6510.02 

673-691 

0.6110.04 

0.5110.01 

0.6410.02 

0.7310.04 

716-725 

0.5210.04 

0.5010.05 

0.6610.12 

770 

0.4510.01 

Mixed  S 

502-543 

0.2710.15 

0.3410.09 

0.2110.07 

559 

0.2710.04 

615-628 

0.2210.04 

0.2910.09 

0.3310.05 

671-689 

0.2410.09 

0.2610.07 

0.3110.03 

0.3610.08 

0.33 

716-726 

0.2610.02 

0.2310.07 

0.2010.05 

0.3010.07 

0.2610.07 

765  0.20±0.10 
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Table  F-2.  Ratios  of  S02  Emission  Index  Values  to  Theoretical  Limits  for  Individual  Engine  Power 
Settings,  Fuel  Sulfur  Loadings  and  Simulated  Altitudes  (Error  Limits  Formed  from 
Standard  Deviations  and  an  Estimate  of  Fractional  Systematic  Error  of  0. 16) 


Power 

Setting, 

T,jc 

Altitude,  m 

16,764 

15,240 

13.716 

12.192 

10.668 

9144 

HighS 

500-538 

0.6310.17 

0.5910.16 

0.7210.17 

572 

0.8310.16 

616-630 

0.8110.17 

0.8410.20 

0.7110.16 

672-686 

0.79±0.19 

0.7310.18 

0.7110.16 

0.9110.20 

0.7310.17 

714-720 

0.65±0.17 

0.6910.18 

0.6810.19 

0.7910.22 

0.7710.17 

783 

0.6710.16 

P-8 

505-506 

0.7210.21 

0.5110.16 

583-619 

0.8010.17 

0.9710.16 

673-691 

0.91±0.17 

0.7510.16 

0.9510.16 

1.0810.17 

716-725 

0.7710.17 

0.7410.18 

0.9810.24 

770 

0.6810.16 

Mixed  S 

502-543 

0.8810.51 

1.1010.34 

0.7010.27 

559 

0.8810.20 

615-628 

0.7110.21 

0.9710.33 

1.0910.24 

671-689 

0.7910.34 

0.8610.28 

1.0310.19 

1.1710.31 

1.1010.16 

716-726 

0.8410.17 

0.7710.28 

0.6310.22 

0.9810.28 

0.8510.29 

765 

0.6610.38 
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Table  F-3.  S02  Emission  Index  Values  Averaged  Over  Engine  Power  Setting,  for  a  Range  of  Fuel 
Sulfur  Loadings  and  Simulated  Altitudes,  with  Standard  Deviations  in 
Averaged  Values  from  Table  F-l. 


Altitude,  m 


Theoretical 

Limit 

16,764 

15,240 

HighS 

2.226 

1.60±0.23 

1.6610.14 

JP-8 

0.672 

0.57±0.06 

Mixed  S 

0.304 

0.25±0.01 

0.2410.02 

13.716 

12.192 

10.668 

9144 

1.5910.20 

1.7110.34 

1.6410.12 

0.5310.02 

0.5410.09 

0.5710.16 

0.2710.05 

0.3310.03 

0.2710.06 

Table  F-4.  Values  of  S02  Fraction  of  Total  Sulfur  in  Exhaust  Averaged  Over  Engine  Power  Setting, 
for  a  Range  of  Fuel  Sulfur  Loadings  and  Simulated  Altitudes,  with  Error  Limits  from  Table  F-3 
Standard  Deviations  and  a  0.16  Estimate  of  Fractional  Systematic  Errors 


Altitude,  m 


Fuel 

16,764 

15,240 

13.716 

12.192 

10.668 

9144 

HighS 

0.7210.19 

0.7410.17 

0.7210.18 

0.7610.23 

0.7410.17 

JP-8 

0.8210.19 

0.7410.16 

0.7810.21 

0.8610.32 

Mixed  S 

0.8210.17 

0.7810.17 

0.8810.23 

1.1010.19 

0.8810.24 

Overall 

Average 

0.7910.18 

0.7610.17 

0.7410.16 

0.8010.23 

0.8810.25 

0.8210.23 
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Table  F-5.  Values  of  S02  Emission  Index  Averaged  Over  Simulated  Altitudes,  for  a  Range  of 
Engine  Power  Settings  and  Fuel  Sulfur  Loadings,  with  Standard  Deviations  in 
Averaged  Values  from  Table  F-l. 


Engine  Power  Setting  Range,  T3  in  °K 

Theoretical 


Limit 

500-545 

560-620 

615-630 

670-690 

710-730 

765-785 

Fuel 

HighS 

2.226 

1.44±0.15 

1.85 

1.7510.15 

1.7210.18 

1.5910.14 

1.49 

JP-8 

0.672 

0.41+0.10 

0.6010.08 

0.6210.09 

0.6310.11 

0.45 

Mixed  S 

0.304 

0.27±0.07 

0.27 

0.2810.06 

0.3010.05 

0.3310.05 

0.20 

Table  F-6.  Values  of  S02  .Fraction  of  Total  Sulfur  in  Exhaust  Averaged  Over  Simulated  Altitudes, 
for  a  Range  of  Engine  Power  Settings  and  Fuel  Sulfur  Loadings,  with  Error  Limits  from  Table  F-3 
Standard  Deviations  and  a  0.16  Estimate  of  Fractional  Systematic  Errors. 


Engine  Power  Setting  Range,  T3  in  °K 


Fuel 

HighS 

442-517 

0.6510.17 

547-654 

0.8310.16 

648-675 

0.7910.16 

JP-8 

0.6210.22 

0.8910.20 

Mixed  S 

0.8910.26 

0.8810.16 

0.9210.25 

Overall 

Average 

0.7310.18 

0.8710.18 

0.8610.22 

749-784 

0.7710.18 

825-850 

0.72+0.17 

917-950 

0.67+0.16 

Overall 

Average 

0.74+0.18 

0.9210.21 

0.8310.21 

0.6810.16 

0.82+0.23 

0.9910.23 

0.8110.20 

0.66+0.16 

0.89+0.23 

0.8910.22 

0.7210.17 

0.6710.16 

0.8210.22 
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Table  F-7.  Bases  for  Systematic  Uncertainty  Estimates 


Estimated  Error  Limit  in  Data  Analysis  Input  Parameter  Fractional  Uncertainty 

in  S02  Emission  Index 

ip 

Exhaust  temperature  ±40  K  from  average  change  in  temperature  from  0.09 

simple  average  to  mixing  fraction  weighted  average 

Exhaust  pressure  ±10  per  cent  from  standard  deviations  in  values  that  should  be  0.05 

the  same  across  exit  plane  (compared  to  ±5  per  cent  from  two  averages  as  above) 

Using  next-nearest  instead  of  nearest  set  of  tuning  coefficients  0.01 

Error  in  baseline,  from  changing  polynomial  order  .  0.05 

Water  line  widths  ±10  per  cent  from  review  article  which  notes  that  0.06 

despite  average  error  bars  from  individual  studies  of  5  per  cent,  comparing 
studies  yields  a  real  uncertainty  of  10  per  cent  or  more 

Sulfur  dioxide  line  widths  ±1 1  per  cent  from  literature  value  of  average  linewidth  0.09 

for  observed  band  of  0.1  ±  0.01 1  wavenumbers  (again,  average  uncertainty  in 
individual  measured  width  values  in  this  source  is  ±0.007  wavenumbers) 

Sulfur  dioxide  band  strength  ±3.2  per  cent,  by  comparing  adjusted  values  from  0.032 

NIST/Orsay  with  HITRAN  values 


The  above,  added  in  quadrature,  result  in  an  estimate  of  fractional  systematic  error  of  0. 16 


A 
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Figure  F-l.  Optical  Layout  of  Tunable  Infrared  Diode  Laser  System. 
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Laser  Intensity,  mV 


Transmission 


Figure  F-3.  Observed  Transmission  Spectrum  (Solid  Line)  and  Single  Temperature  and  Pressure 
Model  Fit  (Dotted  Line). 
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Transmission 


Figure  F-4.  S02  Region  Observed  Transmission  Spectrum  (Solid  Line  and  Circles)  and  Uniform 
Path  (Single  Temperature  and  Pressure)  Model  Fit  By  Data  Analysis  Program  (Dotted  Line).  The  fit 
in  this  spectral  region  is  used  to  determine  the  relative  SOa  column  density. 
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Transmission  -  H20  and  S02  Models  (Dot-Dash  and  Dot)  and  Sum  (Dash) 


1 

0.998 

0.996 

0.994 

0.992 

0.99 

0.988 

0.986 

1332.1  1332.15  1332.2  1332.25  1332.3  1332.35  1332.4  1332.45 

Frequency,  Wavenumbers 

Figure  F-5.  S02  Region  Observed  Transmission  Spectrum  (Solid  Line  and  Circles)  and  Uniform 
Path  (Single  Temperature  and  Pressure)  Model  Spectra  for  the  Individual  Species  S02  (Dotted  Line) 
and  Water  (Dot-Dash  Line)  as  well  as  Sum  of  the  Two  Model  Spectra  (Dashed  Line). 
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